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ABSTRACT 


This  study  investigates  the  self  noise  from  a  ventilated  supercavitating  vehicle.  A 
ventilated  supercavity  is  a  gaseous  envelope  surrounding  an  underwater  vehicle 
that  significantly  reduces  the  drag  felt  by  the  vehicle.  But  the  hydrodynamic  noise 
generated  by  the  creation  of  the  supercavity  could  impact  the  successful  deployment 
of  the  vehicle.  A  principal  source  of  self  noise  for  these  types  of  vehicles  is  sound 
created  by  the  ventilating  gas  jets  impinging  on  the  air-water  interface.  Analytical 
models  of  the  radiated  sound  through  the  interface  have  been  developed.  Sometimes 
jets  impinging  on  the  interface  entrain  bubbles  beneath  the  surface.  This  thesis 
outlines  a  theory  to  predict  the  influence  of  bubbles  near  the  interface. 

Experimental  measurements  were  made  at  the  Naval  Undersea  Warfare  Center 
(NUWC)  in  Newport,  RI  to  test  the  accuracy  of  the  model.  These  measurements 
include  the  unsteady  force  spectrum  of  a  gas  jet  impinging  on  a  rigid  wall.  The 
acoustic  pressure  spectrum  of  a  gas  jet  striking  the  air- water  interface  was  also 
recorded.  The  experimental  results  were  compared  to  theoretical  models  for 
validation. 
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Chapter  1 


Introduction 

A  supercavity  is  a  volume  of  gas  that  completely  encloses  an  underwater  vehicle. 
The  supercavity  is  created  and  sustained  by  continually  injecting  gas  into  the 
region  where  cavitation  first  occurs.  Hydrodynamic  ’self-noise’  generated  by  the 
supercavity  could  interfere  with  the  successful  deployment  of  the  vehicle.  In  order  to 
reduce  this  problem  it  is  necessary  to  understand  the  noise  production  mechanisms. 
One  of  the  principal  sources  of  noise  is  believed  to  be  the  gas  ventilation  system, 
and  the  outcome  of  this  project  will  be  to  study  its  impact  on  the  vehicle  self  noise. 
This  project  is  an  extension  of  work  done  by  Foley  investigating  sound  generated  by 
a  supercavitating  vehicle  [1,2]. 

1.1  Motivation 

The  U.S.  Navy  is  interested  in  developing  high-speed  supercavitating  underwater 
vehicles.  The  Office  of  Naval  Research  is  coordinating  the  collaborative  effort 
between  various  naval  laboratories.  The  primary  goal  of  this  project  is  to 
get  a  better  understanding  of  the  self-noise  generated  by  the  gas  ventilation 
system.  This  includes  investigating  the  influence  of  bubbles  near  the  air-water 
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interface  on  the  radiated  sound.  This  study  will  allow  vehicle  designers  to  make 
informed  decisions  about  the  gas  ventilation  systems  used  in  supercavitating  vehicles. 


1.2  Background 

1.2.1  Cavitation 

Cavitation  is  the  creation  of  pockets  of  air  or  vapor  within  a  liquid  volume. 
Cavitation  often  occurs  when  there  is  a  drop  in  pressure  below  the  liquid’s  vapor 
pressure  causing  it  to  change  phase.  Cavitation  can  occur  in  liquids  that  are  either 
in  motion  or  at  rest.  Pockets  of  cavitation  can  form  both  within  the  body  of  a  liquid 
and  on  solid  boundaries.  Cavitation  on  a  solid  boundary  is  called  fixed  cavitation, 
and  is  of  primary  importance  in  the  study  of  underwater  vehicles.  In  supercavitating 
vehicles,  the  boundary  geometry  of  the  ’cavitator’  at  the  nose  causes  flow  separation 
from  the  solid  wall.  High  velocity  flow  drops  the  pressure,  and  cavitation  occurs 
within  the  separation  region.  These  small  cavities  group  together  to  form  a  larger 
cavity  attached  to  the  boundary.  Without  flow  separation  and  this  region  of  relative 
calm,  any  small  cavities  formed  are  unable  to  attach  to  the  boundary  and  are 
convected  away  by  the  liquid  flow  [3-5]. 

If  a  stable  flow  is  established,  the  cavity  can  either  be  maintained  or  expanded 
downstream  of  the  injection  site.  Cavities  which  begin  and  end  on  the  same  solid 
boundary  are  called  partial  cavities.  Cavities  that  end  downstream  of  the  boundary 
on  which  they  begin  are  called  supercavities.  The  artificial  injection  of  more  gas 
into  an  initial  cavity  in  order  to  expand  or  maintain  it  is  called  ventilation.  The 
role  of  the  gas  ventilation  system  in  the  context  of  supercavitating  vehicles  will  be 
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discussed  in  the  next  section. 


1.2.2  Supercavitating  vehicles 

A  supercavitating  vehicle  is  a  high  speed  underwater  vehicle  that  travels  inside 
a  stable  gaseous  supercavity.  The  elimination  of  friction  with  the  water  reduces 
the  drag  by  up  to  90%  [6,7].  The  supercavity  is  initiated  by  flow  separation  at 
the  cavitator,  which  is  a  specialized  name  for  the  nose  of  these  types  of  vehicles. 
The  cavity  is  then  expanded  by  ventilating  more  gas  into  the  initial  cavity  [3,6,7]. 
The  gas  injection  rate  is  carefully  controlled  to  maintain  cavity  stability  [8,9]. 
Under  stable  high  speed  conditions,  the  gas  exhausts  from  the  rear  of  the  cavity 
py  periodically  shedding  ring  vortices  with  gaseous  cores.  At  lower  speeds,  the 
supercavity  resembles  a  plume  and  gas  escapes  through  two  trailing  hollow  vortex 
tubes  formed  by  buoyancy  induced  bifurcation [3, 10]. 

Supercavitating  vehicles  rely  on  high  frequency  acoustic  sensors  in  the  nose  for 
navigation.  Safe  deployment  of  a  supercavitating  vehicle  therefore  requires  that  any 
interference  from  cavity  self-noise  be  identified,  quantified,  and  possibly  eliminated 
[11-14],  The  principal  sources  of  self-noise  are  thought  to  be  the  gas-jet  ventilation 
system  and  propulsive  sources  at  the  rear  of  the  vehicle  [1,2,15,16].  But  sound 
from  the  propulsive  system  is  confined  primarily  to  lower  frequencies.  Also,  sound 
from  the  rear  of  the  vehicle  is  substantially  attenuated  by  propagation  over  the 
pressure-release  cavity  surface  [17].  Turbulence  in  the  aqueous  boundary  layer  at 
the  edge  of  the  cavitator  generates  sound  and  hydrodynamic  pressure  perturbations 
as  it  convects  across  the  edge  [13,15,18,19].  To  analyze  other  sources  of  sound 
downstream  of  the  cavitator,  it  is  reasonable  to  treat  the  edge  of  the  supercavity 
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as  a  pressure  release  surface.  Then  the  sound  generated  by  turbulence  quadrupoles 
should  contribute  little  to  the  overall  sound,  because  the  pressure  release  surface 
causes  the  pressure  to  vary  as  pwv2M3,  where  the  Mach  number  M  =  v/cw  1,  v 
representing  flow  velocity  and  rhoW)  cw  are  respectively  the  mean  density  and  sound 
speed  in  the  water  [20].  Ventilation  noise  sources  include  direct  gas  impingement  on 
the  cavity  interface  and  the  break-up  of  the  interface  into  a  roughened  spray/bubble 
layer.  These  are  potentially  high  frequency  sources  of  sound  characterized  by 
instabilities  with  length  scales  smaller  than  the  supercavity  radius  [11,12,15,16]. 

Foley  et  al.  [2,6]  used  a  specially  designed  supercavitating  vehicle  to  investigate 
the  radiated  noise  from  cavity  ventilation.  In  a  typical  design  the  ventilating  gas  is 
injected  through  a  series  of  inlet  nozzles  and  redirected  with  a  deflection  device  so 
that  the  gas  strikes  the  interface  at  a  shallow  angle,  producing  a  relatively  smooth 
interface  (Figure  1.1).  In  Foley’s  experiments  (at  the  Penn  State  Applied  Research 
Laboratory  (ARL))  the  deflector  was  removed  so  that  each  jet  impinged  directly 
on  the  interface  at  close  to  normal  incidence.  Impingement  at  normal  incidence  at 
speeds  between  40-50  m/s  leads  to  the  entrainment  of  gas  bubbles  (Figure  1.2). 


5 


Figure  1.1:  Supercavitating  vehicle  in  a  water  tunnel  with  a  smooth  cavity  surface. 


Figure  1.2:  Supercavitating  vehicle  in  a  water  tunnel  with  a  rough  and  bubbly  cavity 
surface. 


It  is  important  to  consider  the  various  paths  of  self  noise  to  the  acoustic  sensors 


at  the  nose  of  the  supercavitating  vehicle.  A  popular  design  choice  is  to  deflect 
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the  gas  jets  to  impact  the  interface  at  a  shallow  angle.  A  potential  issue  with  this 
design  choice  is  that  it  transmits  vibrations  directly  through  the  structure  towards 
the  acoustic  sensors  at  the  nose  (Figure  1.3).  The  total  acoustic  signal  at  the  nose 
would  be  a  combination  of  the  vibrations  transmitted  through  the  structure  and  the 
signal  that  travels  through  the  water.  Supercavitating  vehicles  with  deflectors  tend 
to  produce  a  more  stable  cavity,  which  is  essential  to  the  operation  of  these  vehicles. 
As  a  result,  it  will  be  important  to  understand  the  interaction  between  gas  jets  and 
deflector  plates. 


Figure  1.3:  Acoustic  paths  of  transmission  for  a  supercavitating  vehicle  with  a  de¬ 
flector. 

1.2.3  Free  gas  jets  impinging  on  plates 

A  review  of  existing  literature  showed  that  the  dynamics  of  jet  impingement  are 
well  explored.  Many  experimental  and  theoretical  studies  of  impingement  on  both 
rigid  and  free  surfaces  have  been  made  [21-23].  In  these  studies,  the  heat  transfer 
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aspect  of  the  jet  impingement  was  a  priority  clue  to  its  application  to  manufacturing 
processes.  The  acoustics  of  a  jet  impinging  on  a  plate  have  also  been  explored 
experimentally  [24-26].  These  studies  focus  on  sound  production  mechanisms  like 
shed  vorticity  and  turbulence  rather  than  the  unsteady  motion  of  a  free  surface.  One 
experiment  performed  by  Strong  et  al  [27]  is  more  relevant  to  the  supercavitating 
problem  because  it  measured  the  unsteady  pressure  generated  by  gas  jets  on  a  flat 
surface.  Foley  used  Strong’s  measurements  as  a  check  for  some  of  her  experimental 
measurements  [2],  A  study  of  free  jet  velocities  shows  that  over  small  length  scales, 
the  exit  velocity  of  the  jet  at  the  nozzle  is  the  same  as  the  impact  velocity  of  the  jet 
on  the  interface  [28].  As  a  result,  experiments  in  this  study  were  designed  so  that 
the  jet  exit  velocity  would  be  the  same  as  its  impact  velocity. 


1.2.4  Underwater  Noise 

The  primary  basis  of  this  work  is  Foley’s  previous  investigations  into  sound 
generated  by  the  ventilation  systems  of  supercavitating  vehicles.  Foley  developed  an 
acoustic  model  to  predict  the  radiated  sound  a  jet  impinging  on  a  smooth  air-water 
interface.  This  model  uses  the  measured  unsteady  force  of  a  jet  impinging  on  a 
rigid  wall  as  an  input.  Search  of  the  literature  has  shown  that  there  have  been 
no  underwater  measurements  of  gas  jets  impinging  on  an  air-water  interface.  It 
therefore  became  a  major  goal  of  this  investigation  to  take  underwater  acoustic 
measurements  of  gas  jets  impinging  on  an  air-water  interface.  This  would  provide 
the  first  real  chance  to  compare  the  acoustic  model  to  real  data.  Foley  et  al  [2] 
made  measurements  of  the  unsteady  force  spectrums  of  a  gas  jet  impinging  on  a 
rigid  wall.  These  measurements  will  be  used  as  a  basis  of  comparison  for  the  force 


measurements  described  in  this  thesis. 


The  theories  described  later  in  this  thesis  are  an  extension  of  Foley’s  original 
theory  to  account  for  the  presence  of  bubbles.  This  led  to  considerable  research  into 
the  acoustics  of  bubbly  mixtures.  Wood  predicted  that  the  frequency  dependent 
sound  speed  in  a  bubbly  mixture  is  dependent  on  the  void  fraction  and  resonant 
frequency  of  the  bubbles  in  the  mixture  [29].  The  resonant  frequency  of  bubbles 
was  shown  to  be  a  function  of  the  radius  by  Minnaert[30].  Brennen  provides  an 
excellent  overview  of  the  acoustics  of  bubbly  mixtures[6].  Sound  in  multi-phase 
flows  are  influenced  by  viscous  and  thermal  diffusion  [6,31].  The  results  of  Wood 
and  Minnaert  are  described  in  great  detail,  and  provides  a  good  discussion  of  the 
treatment  of  damping  mechanisms,  called  77,  in  the  frequency  dependent  sound 
speed  equation.  Fox  et  al.  [32]  performed  experiments  to  measure  the  sound  speed 
and  attenuation  through  a  bubbly  fixture.  They  experimentally  determined  that 
a  value  of  77  =  0.5  matched  their  measurements.  But  analytical  predictions  of 
the  attenuation  for  the  experiments  performed  by  Fox  et  al.  [32]  give  a  value  of 
77  =  0.1  [1] .  While  77  is  definitely  dependent  on  the  resonant  frequency  of  the  bubbles, 
it’s  exact  value  is  still  somewhat  ambiguous. 

While  there  have  been  no  underwater  measurements  of  gas  jets  hitting  an  air 
water  interface,  there  have  been  underwater  measurements  of  rainfall [33].  This 
is  not  the  problem  examined  in  this  study,  but  the  measured  results  could  give 
some  insight  into  noise  spectrums  generated  by  impacts  at  the  surface.  Pressure 
spectrums  from  this  study  show  relatively  flat  broadband  signals  from  4-20  kHz. 
Studies  have  also  been  made  of  the  collective  oscillations  of  bubble  clouds  in  salt 
water  formed  by  a  plunging  water  jet  [34],  This  paper  highlights  the  importance  of 
understanding  bubble  dynamics  in  a  saltwater  environment.  Bubbles  are  profuse 
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in  a  saltwater  environment,  and  sound  interacting  with  clouds  of  bubbles  can  be  a 
strong  sound  production  mechanism  [34],  For  random  noise  generated  by  a  gas  jet 
hitting  the  interface,  it  seems  reasonable  to  expect  a  relatively  flat  signal  over  that 
frequency  range.  Analytical  studies  have  also  been  made  to  study  the  impact  noise 
of  a  drop  falling  on  water [35].  This  model  focuses  on  the  impact  of  the  drop  on  the 
interface  as  the  main  mechanism  of  sound  transmission.  The  motivation  for  this 
study  was  to  examine  the  effects  of  spray  from  impinging  gas  jets  on  the  interface  of 
a  supercavity.  These  drops  would  likely  fall  back  down  to  the  interface  and  generate 
noise  without  bubble  entrainment.  Studies  by  Carey  [36]  on  the  acoustics  of  bubble 
plumes  created  by  plunging  water  waves  also  provided  insight  into  the  acoustics 
of  multiphase  fluids.  His  discussion  included  a  summary  of  Wood’s  result,  and 
discussed  the  bubble  cloud  as  an  acoustic  source.  Carey’s  paper  also  described  some 
fundamental  differences  between  bubbles  in  fresh  water  and  salt  water.  Bubbles  in 
salt  water  tend  to  be  smaller  than  those  in  fresh  water  because  the  ionic  properties 
of  the  saltwater  prevent  bubble  coalescence. 


1.3  Approach 

The  first  task  in  preparation  for  taking  measurements  related  to  sound  generated 
by  supercavitating  vehicles  was  to  design  a  surrogate  gas  injection  system.  The 
purpose  of  this  system  was  to  generate  a  steady  jet  of  air  while  recording  the  flow 
data  in  order  to  calculate  the  jet  exit  velocity.  After  creating  a  reliable  system  for 
creating  gas  jets,  this  system  was  ready  to  be  used  in  the  other  experiments  planned 
in  this  study. 

Analytical  predictions  coupled  with  experimental  measurements  have  been 
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made  to  study  the  self  noise  from  ventilated  supercavitating  vehicles.  Foley’s  theory 
of  sound  generated  through  the  air-water  interface  was  extended  to  include  the 
influence  of  a  bubbly  cloud  formation.  The  final  expression  for  radiated  sound 
depends  on  empirical  measurements  of  the  unsteady  force  of  the  gas  jet  impinging 
on  a  rigid  wall.  Force  measurements  were  taken  at  the  Naval  Undersea  Warfare 
Center  (NUWC)  in  Newport,  RI  on  August  14-15,  2012.  Gas  jets  of  several  different 
diameters  were  examined.  The  force  measurements  are  used  as  an  input  to  the 
theoretical  radiated  sound. 

The  sound  of  gas  jets  impinging  on  an  air-water  interface  were  recorded  at  the 
Acoustic  Test  Facility  (ATF)  at  NUWC  on  August  20-21,  2012.  The  measurements 
were  compared  to  Foley’s  original  theory  of  sound  generated  through  the  air-water 
interface.  Some  of  the  gas  jets  were  designed  to  entrain  a  cloud  of  bubbles  beneath 
the  surface.  These  results  were  compared  to  the  predictions  of  sound  generated 
through  a  bubbly  cloud. 


1.4  Outline 

Measurements  of  the  gas  jet  injection  system  are  described  in  Chapter  2.  The 
measurement  of  the  unsteady  force  spectrum  from  gas  jets  impacting  a  rigid  wall  is 
discussed  in  Chapter  3.  The  acoustic  measurements  of  gas  jets  impacting  a  planar 
air-water  interface,  and  a  bubbly  air-water  interface  are  discussed  in  Chapters  4 
and  5.  Chapter  6  describes  in  detail  the  calculation  of  sound  radiated  by  a  gas  jet 
impingement  on  a  gas-water  interface.  Chapter  7  summarizes  the  results  of  the 
previous  chapters  and  provides  suggestion  for  further  research  into  the  self  noise  of 
ventilated  supercavitating  vehicles. 
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Chapter  2 


Measurements  of  the  jet  exit  velocity  and 
in  air  acoustic  spectrum 

This  chapter  discusses  measurements  of  the  gas  jet  injection  system.  Section  2.1  will 
describe  the  experimental  setup  and  procedures.  Section  2.2  will  describe  the  data 
post-processing.  This  will  include  a  description  of  the  fluid  dynamics  theory  used 
to  calculate  the  jet  exit  velocities.  Section  2.3  will  present  the  results  of  the  mea¬ 
surements.  Section  2.4  will  summarize  the  conclusions  based  on  these  measurements. 


2.1  Experimental  setup  and  procedure 

In  order  to  calculate  the  exit  velocity  of  the  gas  jet  at  the  nozzle,  several 
measurements  were  taken  of  the  gas  flowing  through  the  system,  Pressure, 
temperature,  and  flow  rate  sensors  were  used  to  record  p\,  Tj,  and  Q i  respectively 
upstream  of  the  jet  nozzle  (see  Appendix  C  for  sensor  documentation).  A  sketch  of 
the  gas  jet  ventilation  system  is  shown  in  Figure  2.1.  The  nitrogen  gas  is  stored  in 
a  compressed  bottle.  Nitrogen  gas  was  selected  for  this  study  because  it  is  safe  and 
easy  to  obtain.  The  pressure  regulator  controls  the  input  pressure  to  the  nitrogen 
gas,  and  is  used  to  set  the  flow  rate.  The  gas  then  travels  15'  through  a  1/4" 
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high  pressure  air  hose  until  it  reaches  the  sensors.  The  pipe  diameter  through  the 
instrumentation  section  of  the  system  is  1/2".  The  gas  then  travels  another  5'  until 
it  reaches  the  jet  nozzle.  The  instrumentation  was  located  as  close  to  the  jet  nozzle 
as  possible.  Four  jet  nozzles  were  used,  with  exit  diameters  of  1/16",  1/8",  3/16", 
and  0.355".  The  gas  flow  data  was  recorded  using  a  National  Instruments  DAQPad 
attached  to  a  PC  running  LabView  software. 


Temperature 

Sensor 

Pressure 
Sensor 
Flow  Sensor 


Figure  2.1:  Gas  jet  ventilation  system  with  pressure,  temperature,  and  flow  rate 
sensors.  These  sensors  measure  p\ ,  Tj,  and  Q i,  respectively. 


The  gas  jets  were  measured  for  flow  rates  attempting  to  match  the  desired  exit 
velocities  used  in  in  the  acoustic  measurements.  The  goal  was  to  achieve  jet  exit 
velocities  of  220  ft/s,  180  ft/s,  and  140  ft/s.  The  flow  rates  chosen  to  get  close  to 
the  goal  velocities  were  based  on  incompressible  flow  theory.  For  incompressible 
flow,  the  exit  velocity  is  the  flow  rate  divided  by  the  area  of  the  exit  nozzle.  In 
addition  to  this,  measurements  were  made  for  the  1/16"  and  1/8"  jet  nozzles  with 
the  pressure  set  to  50  psi.  These  two  high  velocity  tests  were  meant  to  calculate 
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the  exit  velocity  for  the  jets  that  entrained  a  bubble  cloud  beneath  the  air-water 
interface. 

In  an  effort  to  understand  the  noise  of  the  gas  jet  ventilation  system,  in  air 
recordings  of  the  acoustic  pressure  spectrum  were  taken  of  the  free  jet.  Nozzles  with 
diameters  1/16",  1/8",  3/16",  and  0.355"  were  again  used  for  this  measurement.  Gas 
jets  with  exit  velocities  of  140  ft/s  and  180  ft/s  were  recorded.  The  microphone 
was  positioned  13/8"  from  the  jet  nozzle  (see  Appendix  C  for  microphone  documen¬ 
tation).  This  position  was  chosen  to  match  the  distance  the  gas  jet  would  be  set 
from  the  air-water  interface  in  the  acoustic  experiments.  The  force  spectrum  data 
was  sampled  at  a  rate  of  72  kHz  using  a  National  Instruments  DAQPad  and  a  PC 
running  LabView  software. 


2.2  Data  Analysis 

2.2.1  Compressible  gas  flow  theory 

When  a  fluid  moves  at  speeds  comparable  to  the  speed  of  sound,  the  flow  is 
said  to  be  compressible.  One  of  the  most  important  effects  of  compressibility  is 
choking ,  where  the  flow  rate  is  limited  by  the  speed  of  sound,  c.  Compressible 
flow  calculations  depend  on  both  the  Mach  number  M  =  V/c,  where  V  is  the 
fluid  velocity,  and  the  specific  heat  ratio  7  =  y-  [37].  For  moderate  changes  in 
temperature,  it  is  reasonable  to  use  perfect  gas  assumptions  with  constant  specific 
heats.  This  leads  to  the  perfect  gas  law  p  =  pRT,  where  p  is  pressure,  p  is  density, 
R  is  the  gas  constant,  and  T  is  the  temperature.  The  changes  in  enthalpy,  h,  can  be 
calculated  dh  =  cpdT.  The  gas  injection  system  used  in  ventilated  supercavitating 
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vehicles  can  be  modeled  as  an  adiabatic  isentropic  gas  flow. 

The  mass  flux  through  the  system  is  defined  as  tin  =  peVeAe,  where  pE  is  the 
exit  density,  VE  is  the  exit  velocity,  and  AE  is  the  area  of  the  jet  nozzle. 


VE  = 


m 

PeAe 


(2.2.1) 
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where  p\  is  the  measured  pressure  spectrum  upstream  of  the 
the  exit  pressure  at  the  nozzle.  The  exit  pressure  is  assumed 
atmospheric  pressure.  Substitute  into  (2.2.1) 
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nozzle,  and  pE  is 
to  be  the  same  as 
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From  Bernoulli 
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where  cp  is  the  specific  heat  of  the  gas,  T4  is  the  ambient  temperature,  and  Tj  is  the 
measured  temperature  upstream  of  the  nozzle.  From  the  perfect  gas  law 
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Also  the  mass  flux  fn  =  p\V\A\  =  piQi,  where  Q\  is  the  measured  volumetric  flow 
rate  upstream  of  the  nozzle.  Use  (3.2.5)  to  obtain 
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Substitute  this  expression  into  (2.2.4) 
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1/2 


lApi 


rh(  7  —  1) 


+  2  cvTA 


7^Pi 
777(7  —  1) 


Using  777  =  p\V\A 


(2.2.7) 


—  =  UiA  =  A 

Pi 


7^Pi 
777(7  —  1) 


1 1/2 


+  2cpTa 


jApi 


777(7  —  1) 


(2.2.8) 


Substitute  this  into  (2.2.3) 


n  =  A(n)1/7 

VPs/ 


lApi 

777(7  —  1) 


+  2  c„TA 


1/2 


lApi 


777(7  —  1) 


(2.2.9) 
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This  theory  works  until  the  exit  velocity  of  the  jet  gets  close  to  the  speed  of  sound. 
When  this  happens,  the  flow  speed  is  determined  by  the  critical  pressures  and 
densities  of  the  gas  in  the  system.  The  critical  pressure  is  [37] 


Pc  =  P  1 


7/ (7-1) 


The  critical  density  is  [13] 


Pc 


P  i 


The  velocity  through  the  nozzle  is  then  given  by 


Vc  = 


(2.2.10) 


(2.2.11) 


(2.2.12) 


2.2.2  Post-processing 

The  measurements  pi,  Tj,  and  Q 1  were  used  to  calculate  the  exit  velocity  of  the 
gas  jet  from  the  nozzle.  The  first  step  in  analyzing  the  data  was  adjusting  the  flow 
rate  data  to  match  the  actual  pressure  and  temperature  conditions  in  the  lab.  The 
Sierra  Top-Trak  flow  meter  is  calibrated  to  standard  temperature  ( Tst  =  70°F), 
and  pressure  {psr  =  1  atm)  conditions.  The  formula  for  this  conversion  is 

Qi  =  —  ^Qst  (2.2.13) 

Pi  TSt 

Next  use  equation  (2.2.5)  to  calculate  p\.  Next  calculate  the  mass  flux 


rh  =  piQi 


(2.2.14) 
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After  the  mass  flux,  m,  has  been  calculated,  use  equation  (2.2.7)  to  calculate  the 
velocity  at  the  sensors,  V\.  Finally  plug  V\  into  Equation  (2.2.9)  to  find  the  exit 
velocity  at  the  nozzle,  Ve ■  When  the  exit  velocity  at  the  nozzle  gets  close  to  the 
speed  of  sound,  a  different  set  of  equations  was  used  to  calculate  the  exit  velocity. 
This  only  occurs  in  the  runs  where  the  pressure  at  the  regulator  was  set  to  50  psi. 
In  these  cases,  the  flow  is  choked,  and  the  exit  velocity  is  determined  by  the  critical 
pressures  and  densities  of  the  gas.  Equations  (2.2.10-2.2.12)  are  used  to  calculate 
the  choked  flow  exit  velocity. 

Some  work  needed  to  be  done  on  the  microphone  time  series  data  to  calculate 
the  power  spectral  density.  The  power  spectral  density  of  the  voltage  was  calculated 
for  the  microphone  using  a  Hanning  window  and  N  =  2048  Fourier  transform  points 
using  the  following  formula 


S 


1  2  T 

LFndN 2 


nd 

V  AX* 


r— 1 


(2.2.15) 


where  LF  =  0.375  is  the  Hanning  window  loss  factor,  N  =  2048  is  the  number 
of  Fourier  transform  points,  nd  =  50  is  the  number  of  averages,  X  is  the  Fourier 
transform  of  the  windowed  time  series,  and  X*  is  the  complex  conjugate  of  X.  The 
time  series  used  to  calculate  the  power  spectral  density  began  0.25  s  after  data 
acquisition  began,  and  ended  1.1378  s  later.  A  time  series  of  this  length  ensures 
that  50  averages  are  used  to  calculate  the  power  spectral  density. 


2.2.3  Uncertainty 

The  main  source  of  error  in  the  gas  jet  exit  velocity  is  from  the  flow  rate  sensor  (see 
Appendix  C  for  documentation).  The  Sierra  Top-Trak  flowmeter  has  a  +/—  5% 
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error.  It  was  also  difficult  to  accurately  set  the  pressure  regulator  for  the  low  flow 
rates.  In  hindsight,  a  needle  valve  after  the  pressure  regulator  could  have  been  used 
for  finer  control  of  the  flow  rates. 

The  only  source  of  error  for  the  in-air  measurements  was  the  microphone  itself. 
The  Brucl  and  Kjaer  4938  is  accurate  to  ±1%.  During  processing,  the  raw  voltage 
is  converted  to  decibels.  Due  to  the  logarithmic  nature  of  decibels  and  the  relatively 
small  error  from  the  microphone,  it  is  customary  to  assume  that  the  acoustic 
recordings  are  accurate  to  +/  —  1  dB. 


2.3  Results 

2.3.1  Jet  exit  velocity  measurements 

Plots  of  the  calculated  gas  jet  exit  velocities  for  the  1/16"  nozzle  are  shown  in 
Figures  2. 2-2.4.  It  was  difficult  for  the  gas  jet  injection  system  to  produce  a 
steady  exit  velocity  at  these  low  flow  rates.  This  is  evident  in  the  variability  of  the 
calculated  exit  speed. 
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Jet  Exit  Velocity,  D  =  0.0625",  Goal  V  =  140  ft/s 


Figure  2.2:  Gas  jet  exit  velocity  for  the  1/16"  nozzle  with  a  goal  velocity  of  140  ft/s. 


Jet  Exit  Velocity,  D  =  0.0625",  Goal  V  =  180  ft/s 


Figure  2.3:  Gas  jet  exit  velocity  for  the  1/16"  nozzle  with  a  goal  velocity  of  180  ft/s. 
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Jet  Exit  Velocity,  D  =  0.0625",  Goal  V  =  220  ft/s 


Figure  2.4:  Gas  jet  exit  velocity  for  the  1/16"  nozzle  with  a  goal  velocity  of  220  ft/s. 


Plots  of  the  calculated  gas  jet  exit  velocities  for  the  1/8"  nozzle  are  shown 


in  Figures  2. 5-2. 7.  The  gas  jet  system  was  better  at  controlling  the  higher  flow 
rates.  As  the  nozzle  diameter  and  velocity  increase,  the  calculated  exit  velocity  gets 
smoother. 
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Jet  Exit  Velocity,  D  =  0.125",  Goal  V  =  140  ft/s 


Figure  2.5:  Gas  jet  exit  velocity  for  the  1/8"  nozzle  with  a  goal  velocity  of  140  ft/s. 


Jet  Exit  Velocity,  D  =  0.125",  Goal  V  =  180  ft/s 


Figure  2.6:  Gas  jet  exit  velocity  for  the  1/8"  nozzle  with  a  goal  velocity  of  180  ft/s. 
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Jet  Exit  Velocity,  D  =  0.125",  Goal  V  =  220  ft/s 


Figure  2.7:  Gas  jet  exit  velocity  for  the  1/8"  nozzle  with  a  goal  velocity  of  220  ft/s. 

Plots  of  the  calculated  gas  jet  exit  velocities  for  the  3/16"  nozzle  are  shown  in 
Figures  2.8-2.10.  The  gas  jet  system  maintained  a  relatively  steady  flow  rate  for 
this  nozzle.  The  calculated  exit  velocities  were  slightly  below  the  goal  velocities. 
The  goal  velocities  were  based  on  incompressible  flow  theory.  The  calculated  exit 
velocities  shows  the  difference  between  compressible  and  incompressible  flow  theory 
for  this  case. 
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Jet  Exit  Velocity,  D  =  0.1875",  Goal  V  =  140  ft/s 


Figure  2.8:  Gas  jet  exit  velocity  for  the  3/16"  nozzle  with  a  goal  velocity  of  140  ft/s. 


Jet  Exit  Velocity,  D  =  0.1875",  Goal  V  =  180  ft/s 


Figure  2.9:  Gas  jet  exit  velocity  for  the  3/16"  nozzle  with  a  goal  velocity  of  180  ft/s. 
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Jet  Exit  Velocity,  D  =  0.1875",  Goal  V  =  220  ft/s 


Figure  2.10:  Gas  jet  exit  velocity  for  the  3/16"  nozzle  with  a  goal  velocity  of  220 
ft/s. 


Plots  of  the  calculated  gas  jet  exit  velocities  for  the  0.35"  nozzle  are  shown  in 
Figures  2.11-2.13.  The  gas  jet  system  was  able  to  maintain  a  steady  flow  rate  for 
this  nozzle.  The  calculated  exit  velocities  were  much  lower  than  the  goal  velocities. 
The  results  from  this  nozzle  indicate  that  there  is  a  greater  difference  between 
compressible  and  incompressible  flow  theory  when  the  nozzle  size  is  increased. 
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Jet  Exit  Velocity,  D  =  0.355",  Goal  V  =  140  ft/s 


KL 


10  20  30  40  50 

Time,  s 


Figure  2.11:  Gas  jet  exit  velocity  for  the  0.355"  nozzle  with  a  goal  velocity  of  140 
ft/s. 


Jet  Exit  Velocity,  D  =  0.355",  Goal  V  =  180  ft/s 


Figure  2.12:  Gas  jet  exit  velocity  for  the  0.355"  nozzle  with  a  goal  velocity  of  180 
ft/s. 
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Jet  Exit  Velocity,  D  =  0.355",  Goal  V  =  220  ft/s 


Figure  2.13:  Gas  jet  exit  velocity  for  the  0.355"  nozzle  with  a  goal  velocity  of  220 
ft/s. 


When  the  pressure  regulator  was  set  to  50  psi  for  the  1/16"  and  1/8"  nozzles, 
the  jet  exit  velocity  gets  close  to  the  speed  of  sound.  The  critical  exit  velocity  is 
calculated  using  Equation  2.2.12.  Plots  of  the  critical  velocity  for  the  1/16"  nozzle 
are  shown  in  Figure  3.14.  For  both  nozzles  the  critical  velocity  was  ~  1056  ft/s. 
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Figure  2.14:  Gas  jet  exit  velocity  for  the  1/16"  nozzle  with  the  pressure  regulator 
set  to  50  psi. 


2.3.2  In  air  acoustic  pressure  spectrum  measurements 

The  acoustic  power  spectrum  for  the  1/16"  nozzle  with  an  exit  velocity  of  180 
ft/s  is  shown  in  Figure  2.15.  The  biggest  spike  in  the  spectrum  occurs  near  7 
kHz.  Some  smaller  peaks  occur  at  around  16  and  24  kHz.  The  acoustic  power 
spectrum  for  the  1/8"  nozzle  with  an  exit  velocity  of  180  ft/s  is  plotted  in  Figure 
2.16.  This  spectrum  also  shows  a  spike  near  7  kHz,  but  is  otherwise  mostly 
flat.  The  acoustic  power  spectrum  for  the  3/16"  nozzle  with  an  exit  velocity  of 
180  ft/s  is  shown  in  Figure  2.17.  This  spectrum  shows  the  same  spike  near  7 
kHz.  The  acoustic  power  spectrum  for  the  0.355"  nozzle  with  an  exit  velocity  of 
180  ft/s  is  shown  in  Figure  2.18.  This  spectrum  appears  quite  different  from  the 
other  nozzles.  It  has  quite  a  few  more  spikes,  including  some  large  ones  above  20  kHz. 
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Microphone  Power  Spectral  Density,  D  =  1/16",  V  =  180  ft/s 


Figure  2.15:  In  air  measurement  of  the  free  gas  jet  with  a  nozzle  diameter  of  1/16" 
and  an  exit  velocity  of  180  ft/s. 


Microphone  Power  Spectral  Density,  D  =  1/8",  V  =  180  ft/s 


Figure  2.16:  In  air  measurement  of  the  free  gas  jet  with  a  nozzle  diameter  of  1/8" 
and  an  exit  velocity  of  180  ft/s. 
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Microphone  Power  Spectral  Density,  D  =  3/16",  V  =  180  ft/s 


Figure  2.17:  In  air  measurement  of  the  free  gas  jet  with  a  nozzle  diameter  of  3/16" 
and  an  exit  velocity  of  180  ft/s. 


Microphone  Power  Spectral  Density,  D  =  0.355",  V  =  180  ft/s 


Figure  2.18:  In  air  measurement  of  the  free  gas  jet  with  a  nozzle  diameter  of  0.355" 
and  an  exit  velocity  of  180  ft/s. 


31 


2.4  Conclusions 

The  measurements  of  the  gas  jet  exit  velocities  are  important  in  the  analysis  of  the 
acoustic  data.  These  velocities  will  be  used  to  try  and  determine  the  relationship 
between  jet  exit  velocity  and  radiated  acoustic  sound.  The  experimental  measure¬ 
ments  of  the  gas  jet  exit  velocities  from  the  gas  injection  system  show  that  the  exit 
velocity  is  not  always  steady.  This  is  especially  true  at  the  lower  flow  rates  used  in 
the  smaller  nozzles. 

The  in  air  recordings  of  the  acoustic  spectrum  of  the  free  gas  jet  show  that 
there  are  some  tones  generated  by  the  jet  injection  system.  Spectrum  for  the 
1/16",  1/8",  and  3/16"  nozzles  all  show  a  spike  at  7  kHz.  This  tone  appears  to  be  a 
characteristic  feature  of  the  gas  jet  injection  system  for  these  nozzles.  This  could  be 
a  result  of  irregularities  in  the  flow  through  the  system.  An  effort  was  made  to  keep 
the  flow  as  smooth  as  possible  within  the  system,  but  the  attachment  of  the  sensors 
to  measure  the  flow  rate  created  changes  in  the  interior  pipe  diameter. 


Chapter  3 


Measurements  of  the  unsteady  force 
spectrum  of  a  gas  jet  impacting  a  rigid 
wall 

This  chapter  describes  the  force  spectrum  measurements  performed  at  NUWC 
Newport,  RI  on  August  14-15,  2012.  Section  3.1  will  describe  the  experimental 
setup,  including  the  sensors  used  and  data  acquisition  parameters.  It  also  will 
outline  the  procedure  used  to  collect  the  force  spectrum  data.  Section  3.2  will 
describe  the  data  post-processing,  as  well  as  the  uncertainty  of  the  measured  data. 
Section  3.3  will  present  the  recorded  measurements,  while  section  3.4  will  compare 
these  results  to  the  measurements  previously  made  by  Foley.  Section  3.5  will 
highlight  the  important  conclusions  from  this  experiment. 


3.1  Experimental  setup  and  procedure 

The  force  measurement  experiments  took  place  in  an  anechoic  chamber  located  at 
NUWC.  The  chamber  is  designed  to  have  as  little  background  noise  as  possible, 
and  is  completely  lined  with  sound  dampening  tiles  to  virtually  eliminate  reflected 
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noise.  An  array  of  force  sensors  was  mounted  to  a  massive  rigid  plate  to  record 
the  unsteady  force  spectrum  of  the  gas  jet  hitting  a  rigid  wall.  An  array  of  7  PCB 
Piezotronics  force  sensors  was  attached  to  12"  square  steel  plates  as  shown  in  Figures 
3. 1-3. 3.  A  PCB  Piezotronics  accelerometer  (see  Appendix  C  for  documentation)  was 
attached  to  the  steel  plates  to  measure  the  vibrations  of  the  entire  structure.  These 
steel  plates  were  mounted  to  a  translation  table  that  could  accurately  position 
the  rigid  wall  near  the  gas  jet.  The  data  from  6  of  the  force  sensors  and  the 
accelerometer  was  acquired  using  a  National  Instruments  DAQPad  using  Labview 
software  at  a  sampling  rate  of  90  kHz.  The  force  data  was  passed  through  a  low 
pass  filter  set  at  40  kHz,  and  then  20  dB  of  gain  was  applied  to  the  signal. 


Layout  of  Force  Sensors 


Figure  3.1:  Layout  of  the  force  sensors  on  the  steel  plate. 
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Figure  3.2:  Top  view  of  the  force  sensors  and  accelerometer. 


Figure  3.3:  Side  view  of  the  force  sensors  and  jet  nozzle. 


The  model  of  the  gas  injection  system  investigated  in  these  tests  is  a  single  jet 
of  gas  hitting  a  rigid  wall  at  a  normal  angle.  In  order  to  create  these  conditions, 
a  system  was  devised  that  could  shoot  the  air  jets  at  controllable  and  measurable 
velocities.  The  gas  used  in  the  tests  came  from  a  compressed  Nitrogen  bottle. 
Nitrogen  was  chosen  for  these  experiments  because  of  its  stability.  The  flow  rate 
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was  controlled  with  a  0-250  psi  pressure  regulator.  The  flow  rate  was  measured 
using  a  Sierra  Top-Trak  flowmeter  (see  Appendix  C  for  documentation).  A  National 
Instruments  DAQPad  attached  to  a  PC  running  Labview  software  acquired  the 
data  from  the  flowmeter. 

After  passing  through  the  flowmeter,  the  gas  flowed  through  a  1/4"  air  hose  to 
detachable  nozzles  positioned  next  to  the  rigid  wall.  The  diameters  of  the  nozzles 
used  in  the  tests  were  1/16",  1/8",  3/16",  and  0.355",  and  the  nominal  exit  velocities 
investigated  were  220  ft/s,  180  ft/s,  and  140  ft/s.  There  were  two  additional  runs 
using  the  1/16"  and  1/8"  nozzles  with  the  input  pressure  set  to  50  psi.  These  were 
chosen  to  gather  force  spectrum  measurements  for  the  bubble  entraining  jets.  The 
jet  was  positioned  1.375"  from  the  force  sensor  array.  Table  3.1  summarizes  the  test 
parameters. 


Nozzle 

Diameter 

Gas  Jet  Exit 
Velocities 

Pressure  Settings  for 
Bubble  Entrainment 

Jet  Distance 
from  Interface 

1/16” 

140,  180,  220  fps 

50  psi 

1  3/8” 

1/8” 

140,  180,  220  fps 

50  psi 

1  3/8” 

3/16” 

140,  180,  220  fps 

1  3/8” 

0.355” 

140,  180,  220  fps 

1  3/8” 

Table  3.1:  Test  matrix  for  the  force  measurement  experiments. 


The  procedure  for  acquiring  the  test  data  was  the  same  in  each  trial. 
First,  the  pressure  regulator  was  adjusted  until  the  flow  rate  meter  displayed 
the  appropriate  reading  to  achieve  the  desired  nominal  exit  velocity.  This  flow 
rate  was  calculated  determine  by  an  incompressible  flow  approximation  and  the 
nozzle  exit  diameter.  For  the  cases  where  the  input  pressure  was  set  to  50  psi, 
a  shutoff  valve  in  front  of  the  flowmeter  was  closed.  The  regulator  then  set 
the  system  pressure  to  50  psi,  and  the  valve  was  reopened.  After  the  gas  was 
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flowing  at  the  desired  rate,  the  Labview  software  was  used  to  begin  the  data  acqui¬ 
sition.  After  5  seconds,  the  data  acquisition  was  stopped,  and  the  gas  flow  turned  off. 


3.2  Data  Analysis 


3.2.1  Post  processing 


The  power  spectral  density  of  the  voltage  was  calculated  for  the  force  sensors  and 
accelerometer  using  a  Hanning  window  and  N  =  2048  Fourier  transform  points 
using  the  following  formula 


S 


1  2  T 

LFndN 2 


nd 

V  xx* 


r=  1 


(3.2.1) 


where  LF  =  0.375  is  the  Hanning  window  loss  factor,  N  =  2048  is  the  number 
of  Fourier  transform  points,  nd  =  50  is  the  number  of  averages,  X  is  the  Fourier 
transform  of  the  windowed  time  series,  and  A"*  is  the  complex  conjugate  of  X.  The 
time  series  used  to  calculate  the  power  spectral  density  began  0.25  s  after  data 
acquisition  began,  and  ended  1.1378  s  later.  A  time  series  of  this  length  ensures 
that  50  averages  are  used  to  calculate  the  power  spectral  density. 


3.2.2  Uncertainty 

The  main  source  of  inaccuracy  in  this  measurement  was  the  vibration  of  the 
entire  plate  structure.  An  accelerometer  was  attached  to  the  plate  to  measure  the 
vibrations.  Another  factor  to  consider  is  the  potential  changes  in  velocity  inherent 
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in  the  gas  jet  system.  The  flow  measurements  indicate  that  the  exit  velocity  is 
steadier  at  the  higher  flow  rates.  Finally,  the  force  sensors  themselves  have  a  +/- 
1%  error  (see  Appendix  C  for  documentation). 

3.3  Results 

The  power  spectral  densities  for  the  1/16"  diameter  nozzle  are  shown  in  Figure  3.4. 
It  was  difficult  to  accurately  adjust  the  pressure  regulator  to  achieve  the  desired  exit 
velocities  for  the  1/16"  nozzle.  For  the  1/16"  nozzle,  the  data  hits  the  noise  floor  at 
8  kHz.  Figures  3. 5-3. 6  show  the  power  spectral  densities  measured  from  the  center 
sensor  for  the  1/8"  and  3/16"  nozzles  respectively.  The  power  spectral  densities  dip 
between  1-4  kHz,  rise  back  up  3  dB,  and  then  descend.  Figure  3.7  shows  the  power 
spectral  density  measured  from  the  center  sensor  for  the  .355"  nozzle.  The  power 
spectral  densities  for  the  .355"  nozzle  do  not  have  a  local  minimum  between  1-4 
kHz.  Consistent  features  in  the  data  are  spikes  at  15  kHz  and  32  kHz.  The  spike  at 
32  kHz  can  be  attributed  to  the  natural  resonance  frequency  of  the  force  sensors 
(see  Appendix  C  for  sensor  documentation). 

Figure  3.8  plots  the  power  spectral  densities  measured  from  all  the  force  sensors 
as  well  as  the  accelerometer  data  for  the  1/8"  nozzle  at  222  ft/s.  The  graph  on  the  left 
is  the  force  measurements.  The  blue  curve  corresponds  to  the  All  force  sensor  which 
was  mounted  in  the  center.  The  other  curves  represent  the  outer  ring  of  sensors.  The 
graph  on  the  right  is  the  power  spectral  density  calculated  from  the  accelerometer. 
A  spike  in  the  accelerometer  data  at  15  kHz  indicates  the  entire  structure  was 
vibrating  at  this  frequency.  This  could  be  related  to  the  spike  in  force  data  at  15  kHz. 


Figure  3.9  plots  the  power  spectral  densities  measured  from  the  full  array  of 
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force  sensors  and  the  accelerometer  for  the  1/16"  nozzle  with  pressure  set  to  50 
psi.  The  force  power  spectral  density  is  relatively  flat  for  the  All  sensor,  with  the 
exception  of  the  resonant  peak  at  32  kHz.  The  outer  force  sensors  have  a  spike 
between  7-8  kHz.  The  accelerometer  shows  that  the  entire  structure  is  vibrating 
at  frequencies  greater  than  7  kHz.  Figure  3.10  plots  the  power  spectral  densities 
measured  from  the  full  array  of  force  sensors  and  the  accelerometer  for  the  1/8" 
nozzle  with  the  pressure  set  to  50  psi.  The  force  power  spectral  density  is  flat  for 
the  All  sensor  until  10  kHz,  where  it  decreases  by  6  dB  and  then  flattens  out  again. 
The  outer  force  sensors  have  a  spike  between  7-10  kHz.  The  accelerometer  data  for 
the  1/8"  nozzle  is  similar  to  the  data  from  the  1/16"  nozzle  with  slightly  higher 
amplitudes. 


Force  Sensor  Power  Spectral  Density,  D  =  1/16" 


Figure  3.4:  Power  spectral  density  measurements  of  the  force  sensors  for  1/16"  jet 
diameter. 
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Force  Sensor  Power  Spectral  Density,  D  =  1/8" 


Figure  3.5:  Power  spectral  density  measurements  of  the  force  sensors  for  1/8"  jet 
diameter. 


Force  Sensor  Power  Spectral  Density,  D  =  3/16" 


Figure  3.6:  Power  spectral  density  measurements  of  the  force  sensors  for  3/16"  jet 
diameter. 
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Force  Sensor  Power  Spectral  Density,  D  =  .355" 


Figure  3.7:  Power  spectral  density  measurements  of  the  force  sensors  for  0.355"  jet 
diameter. 


Figure  3.8:  Force  sensor  power  spectral  density  for  the  1/8"  diameter  nozzle  and  exit 
velocity  of  222  ft/s.  The  plot  to  the  right  is  the  accelerometer  data  from  the  same 


run. 
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Figure  3.9:  Force  sensor  power  spectral  density  for  the  1/16"  diameter  nozzle  with 
the  input  pressure  set  to  50  psi.  The  plot  to  the  right  is  the  accelerometer  data  from 
the  same  run. 


Figure  3.10:  Force  sensor  power  spectral  density  for  the  1/8"  diameter  nozzle  with 
the  input  pressure  set  to  50  psi.  The  plot  to  the  right  is  the  accelerometer  data  from 
the  same  run. 
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3.3.1  Comparison  with  Foley 

Foley  et  al  performed  a  similar  set  of  measurements  of  the  unsteady  force  of  gas 
jets  impinging  on  a  rigid  plate  [1,2].  The  difference  in  the  two  experiments  was 
the  gas  jet  injection  system.  Foley  used  a  modified  test  vehicle  to  generate  the  gas 
jets[2].  This  turned  out  to  be  a  significant  source  of  vibration  that  interfered  with 
the  measurement.  The  gas  jet  system  used  in  this  experiment  used  a  system  of  high 
pressure  hoses  that  could  deliver  a  measurable  amount  of  gas  without  unnecessary 
vibrations.  Figure  3.11  plots  a  comparison  of  the  measured  unsteady  force  spectrum 
of  the  two  experiments  for  a  1/8"  nozzle  with  an  exit  velocity  of  140  ft/s.  Foley’s 
data  appears  to  be  noisier  between  3-10  kHz,  and  the  measurements  from  this 
experiment  show  more  of  the  sensor  resonant  spike  at  32  kHz.  Other  than  that  the 
two  measurements  are  similar.  Figure  3.12  shows  shows  the  comparison  of  the  force 
spectrums  for  a  1/8"  nozzle  with  an  exit  velocity  of  180  ft/s.  Foley’s  data  appears 
to  be  much  noisier  between  3-20  kHz.  But  just  before  Foley’s  frequency  jumps  up, 
it  lines  up  well  with  the  measurements  made  in  this  experiment. 
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Comparison  of  Force  Pressure  Spectrums  Made  By  Foley  and  Gilbert,  140  ft/s 


Figure  3.11:  Comparison  with  the  unsteady  force  measurements  made  by  Foley  for 
an  1/8"  nozzle  with  an  exit  velocity  of  140  ft/s. 


Comparison  of  Force  Pressure  Spectrums  Made  By  Foley  and  Gilbert,  180  ft/s 


Figure  3.12:  Comparison  with  the  unsteady  force  measurements  made  by  Foley  for 
an  1/8"  nozzle  with  an  exit  velocity  of  180  ft/s. 


Foley  [1]  estimated  the  form  of  from  measurements  of  the  surface  force 

spectrum  produced  by  low  speed  jet  impact  on  a  nominally  rigid  wall.  From  these 
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measurements  it  was  deduced  that  &ff(w)  is  well  approximated  by  the  empirical 
formula 


F^D 


a 


+  e2(uD/U) 


(3.3.1) 


where  a  and  e  ~  1.6  are  constants.  The  value  of  a  is  adjusted  to  give  the  best 
fit  to  the  measurements,  and  assumes  a  mean  value  of  about  20  for  Foley’s  [1] 
measurements.  Fa  is  the  nominal  mean  force  exerted  by  the  jet  on  the  wall, 
and  F2  =  (p^f/2)2(|Zl2)2  where  pa  is  the  mean  gas  density  (~  1.2 kg/m3  in  the 
experiments) . 


This  empirical  formula  was  compared  to  the  measured  force  spectrums  from  this 
experiment  for  the  1/8"  nozzle  in  order  to  be  consistent  with  Foley’s  experiments. 
Figure  3.13  compares  the  normalized  force  pressure  spectrum  with  the  empirical 
curve  fit  for  the  1/8"  nozzle  at  an  exit  velocity  of  220  ft/s.  In  this  case  the  value  of  a 
was  60  and  e  was  0.7.  Figure  3.14  compares  the  normalized  force  pressure  spectrum 
with  the  empirical  curve  fit  for  the  1/8"  nozzle  at  an  exit  velocity  of  180  ft/s.  In 
this  case  the  value  of  a  was  60  and  e  was  0.7.  Figure  3.15  compares  the  normalized 
pressure  spectrum  with  the  empirical  curve  fit  for  the  1/8"  nozzle  at  an  exit  velocity 
of  140  ft/s.  In  this  case  the  value  of  a  was  100  and  e  was  0.5.  The  changing 
values  of  a  and  e  indicate  that  these  parameters  are  dependent  on  the  jet  exit  velocity. 
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Force  Measurement  Empirical  Curve  Fit  for  D  =  1/8"  at  220  ft/s 


Figure  3.13:  Comparison  of  the  normalized  unsteady  force  for  the  1/8"  nozzle  at  an 
exit  velocity  of  220  ft/s  with  the  empirical  curve  fit  from  equation  3.3.1,  a  =  60,  e 


=  0.7 


Force  Measurement  Empirical  Curve  Fit  for  D  =  1/8"  at  180  ft/s 


Figure  3.14:  Comparison  of  the  normalized  unsteady  force  for  the  1/8"  nozzle  at  an 
exit  velocity  of  180  ft/s  with  the  empirical  curve  fit  from  equation  3.3.1,  a  =  60,  e 


=  0.7 
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Force  Measurement  Empirical  Curve  Fit  for  D  =  1/8"  at  140  ft/s 


Figure  3.15:  Comparison  of  the  normalized  unsteady  force  for  the  1/8"  nozzle  at  an 
exit  velocity  of  140  ft/s  with  the  empirical  curve  fit  from  equation  3.3.1,  a  =  100,  e 
=  0.7 


3.4  Conclusions 

The  analytical  model  for  radiated  acoustic  noise  depends  on  using  an  experimentally 
measured  unsteady  force  spectrum  as  an  input.  Great  care  was  made  during 
the  design  of  this  experiment  to  eliminate  unwanted  noise  due  to  vibrations. 
Comparisons  with  Foley’s  data  show  generally  good  agreement  with  less  vibrational 
noise.  The  consistency  between  the  two  measurements  is  encouraging,  and  indicates 
that  the  unsteady  force  spectrum  measured  in  this  experiment  should  serve  as  a 
good  input  to  the  acoustic  model. 


Chapter  4 


Acoustic  measurements  of  gas  jets 
impacting  a  planar  air-water  interface 

This  chapter  involves  the  acoustic  measurements  of  the  underwater  noise  from  gas 
jets  impinging  on  the  air- water  interface.  Section  4.1  describes  the  experimental 
setup  and  procedure  for  taking  the  measurements.  Section  4.2  contains  descriptions 
of  the  data  analysis,  including  any  post-processing,  and  a  discussion  of  sources  of 
uncertainty  in  the  data.  Section  4.3  presents  the  results  in  graphical  format  and 
discusses  some  of  the  trends  in  the  data.  Section  4.4  summarizes  the  important 
conclusions  drawn  from  the  experiments. 


4.1  Experimental  setup  and  procedures 

Experiments  were  performed  at  the  Acoustic  Test  Facility  at  NUWC  to  measure  the 
noise  generated  by  gas  jets  hitting  the  air-water  interface  Aug.  20-21  of  2012  at  the 
Acoustic  Test  Facility  (ATF)  in  NUWC  Newport,  RI.  The  ATF  is  a  large  indoor 
laboratory  that  specializes  in  accurate  underwater  sound  recordings.  The  facility 
features  a  700,000  gallon  test  tank,  with  dimensions  of  60  ft.  long,  40  ft.  wide,  and 
35  ft.  deep.  Portions  of  the  tank  are  lined  with  acoustic  wedges  to  dampen  out 
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reflections. 

The  gas  ventilation  system  was  positioned  on  a  catwalk  above  the  acoustic 
tank.  The  jet  nozzle  was  attached  to  a  lowering  arm  and  positioned  the  jet  1.75" 
above  the  interface.  A  total  of  4  hydrophones  were  positioned  in  the  tank  beneath 
the  catwalk  to  record  the  underwater  sound  as  shown  in  Figure  4.1.  The  primary 
hydrophone  was  oriented  horizontally  and  positioned  3  m  directly  below  the  gas 
jet.  Another  2  hydrophones  were  oriented  vertically  and  positioned  1  m  away  from 
the  jet  and  1  nr  below  the  surface.  The  fourth  hydrophone  was  again  oriented 
horizontally  and  positioned  1  m  away  from  the  jet  and  1  m  below  the  surface. 
The  orientation  of  the  hydrophones  is  important  because  of  the  directivity  of  the 
hydrophones.  The  hydrophones  are  least  sensitive  in  the  direction  towards  the  cable 
they  are  attached  to.  The  fourth  hydrophone  was  added  so  there  would  be  two 
phones  oriented  in  the  same  direction  for  comparison. 
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Cal  walk 


H3: 1  m  Offset,  1  m 
Depth.  Vertical 
Orientation 


Gas  Jet 


fJl  HI:  0  m  Offset.  3  m 
“  Down,  Horizontal 
Orientation 


_  H2: 1  m  Offset,  1  m 

Depth,  Vertical 
Orientation 


H4: 1  m  Offset.  1  m 
Down,  Horizontal 
Orientation 


Figure  4.1:  Hydrophone  layout  for  acoustic  experiments  on  Aug.  20-21  of  2012  at 
the  ATF. 

Data  was  acquired  from  the  hydrophones  with  the  data,  acquisition  system 
at  the  ATF.  The  acoustic  data  was  sampled  at  a  rate  of  150  kHz.  The  time 
series  was  then  converted  to  a  power  spectrum  using  16,384  points,  and  the 
final  spectrum  recorded  was  a  result  of  16  spectral  averages.  In  this  experi¬ 
ment,  several  jet  diameters  and  jet  velocities  were  examined.  Jet  diameters 
of  1/16",  1/8",  3/16"  and  0.355"  at  exit  velocities  of  140  ft/s,  180  ft/s,  and  220 
ft/s  were  investigated.  Table  4.2  summarizes  the  test  matrix  used  in  this  experiment. 


For  each  trial,  the  following  procedure  was  used  to  acquire  the  data.  First, 
recording  of  the  gas  instrumentation  sensors  was  started  in  Labview.  Then  the  gas 
flow  was  turned  on,  and  the  pressure  regulator  was  adjusted  to  get  the  gas  to  the 
desired  flow  rate.  This  flow  rate  was  estimated  using  incompressible  flow  theory, 
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Nozzle 

Diameter 

Gas  Jet  Exit 
Velocities 

Pressure  Settings  for 
Bubble  Entrainment 

Jet  Distance 
from  Interface 

1/16” 

140,  180,  220  fps 

50  psi 

1  3/8” 

1/8” 

140,  180,  220  fps 

50  psi 

1  3/8” 

3/16” 

140,  180,  220  fps 

1  3/8” 

0.355” 

140,  180,  220  fps 

1  3/8” 

Table  4.1:  Test  matrix  for  acoustic  measurements  on  Aug.  20-21  of  2012  at  ATF. 


along  with  the  nozzle  diameter  and  desired  exit  velocity.  After  the  desired  flow  rate 
was  achieved,  the  acoustic  data  acquisition  system  was  started.  After  the  system 
had  performed  the  16  spectral  averages,  gas  flow  was  shut  off,  followed  by  the  gas 
data  acquisition  system. 


4.2  Data  analysis 

4.2.1  Post-processing 

Data  was  acquired  from  the  hydrophones  with  the  data  acquisition  system  at 
the  ATF.  The  acoustic  data  was  sampled  at  a  rate  of  150  kHz.  The  time  series 
was  then  converted  to  a  power  spectrum  using  16,384  points,  and  the  final 
spectrum  recorded  was  a  result  of  16  spectral  averages.  The  data  was  then 
smoothed  over  frequency  using  a  running  average  with  a  frequency  window  of  82  Hz 
and  90  percent  overlap.  Figure  4.2  shows  the  results  of  applying  this  filter  to  the  data. 
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Figure  4.2:  Comparisons  of  smoothed  pressure  frequency  spectrum  and  raw  pressure 
frequency  spectrum  data. 

4.2.2  Uncertainty 

The  principal  sources  of  uncertainty  in  this  experiment  are  from  the  gas  jet  injection 
system.  The  measurements  of  the  gas  jet  exit  velocities  shows  that  the  exit 
velocity  can  fluctuate,  especially  at  lower  flow  rates.  This  could  impact  the  quality 
of  data  for  the  lower  flow  rates,  especially  the  gas  jets  using  the  1/16"  nozzle. 
Measurements  of  the  in  air  pressure  spectrum  of  the  free  jet  shows  that  the  system 
has  a  characteristic  tone  at  7  kHz  for  the  1/16",  1/8",  and  3/16"  nozzles.  Individual 
bubbles  could  also  potentially  be  contributing  to  the  acoustic  spectrum.  Bubbles 
on  the  order  of  100  microns  could  be  present  near  the  interface,  undetectable  by 
the  naked  eye.  The  only  other  source  of  error  could  be  from  the  hydrophones.  The 
acoustic  data  is  converted  to  decibels,  and  due  to  their  logarithmic  nature  the  error 
is  customarily  marked  as  +/  —  1  dB. 
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4.3  Results 

Figures  4. 3-4. 6  show  the  smoothed  sound  pressure  spectrum  for  all  of  the  trials  with 
no  visible  bubble  entrainment.  The  data  in  Figure  4.3  shows  very  little  variation  in 
the  sound  pressure  levels  as  a  function  of  velocity  for  the  1/16"  nozzle.  The  data 
in  Figure  4.4  shows  a  dependence  of  the  sound  pressure  level  as  a  function  of  the 
velocity  for  the  1/8"  jet  nozzle.  There  also  appear  to  be  relative  maximas  near  16 
and  23  kHz  for  this  case.  The  data  in  Figure  4.5  shows  variation  between  the  runs 
at  different  velocities  for  the  3/16"  jet.  The  spectrum  also  becomes  pretty  flat  at 
the  higher  frequencies  for  this  nozzle.  The  data  in  Figure  4.6  shows  quite  a  bit  of 
variation  between  the  velocities  for  the  1/4"  jet.  For  velocities  above  168  ft/s  a  peak 
develops  between  40-45  kHz.  But  when  the  jet  velocity  is  reduced  to  158  ft/s  the 
peak  vanishes. 


Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/16" 


Figure  4.3:  Measured  acoustic  pressure  spectrums  for  the  1/16"  nozzle  jets. 
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Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/8" 
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Figure  4.4:  Measured  acoustic  pressure  spectrums  for  the  1/8"  nozzle  jets. 


Hydrophone  Pressure  Power  Spectral  Density,  D  =  3/16" 


Figure  4.5:  Measured  acoustic  pressure  spectrums  for  the  3/16"  nozzle  jets. 
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Hydrophone  Pressure  Power  Spectral  Density,  D  =  0.355" 


Figure  4.6:  Measured  acoustic  pressure  spectrums  for  the  0.355"  nozzle  jets. 

Pictures  were  taken  of  gas  jets  impinging  on  the  air-water  interface  in  a  small 
water  tank.  This  tank  would  serve  well  for  imaging  purposes,  but  did  not  serve 
well  as  an  acoustic  tank.  Pictures  were  taken  of  gas  jets  impinging  on  the  air-water 
interface  with  nozzles  of  1/16"  and  1/8"  with  exit  velocities  of  140,  180,  and  220 
ft/s.  The  images  were  taken  with  a  high  speed  camera  sampling  at  a  rate  of  500  Hz. 

Figure  4.7  shows  the  indentation  created  by  a  1/16"  nozzle  jet  with  an  exit 
velocity  of  220  ft/s.  In  this  picture  the  indentation  is  at  a  relative  maximum,  and 
penetrates  about  3/8"  below  the  surface.  But  the  indentation  is  not  fixed  in  time. 
Figure  4.8  shows  another  image  from  the  same  recording  5  frames  later.  In  this 
picture  the  indentation  has  decreased  in  size.  Figure  4.9  shows  the  indentation 
created  by  a  1/8"  nozzle  at  an  exit  velocity  of  220  ft/s.  The  indentation  is  at  a 
maximum  of  about  1.25"  in  this  picture.  There  are  also  several  bubbles  visible  off 
to  the  side  of  the  jet  impact  area.  These  are  likely  from  water  splashing  back  down 
onto  the  interface.  The  impact  of  the  jet  on  the  water  surface  sent  some  droplets 
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up  into  the  air  which  then  fell  back  onto  the  interface.  Figure  4.10  shows  another 
image  of  the  1/8"  nozzle  at  220  ft/s  8  frames  after  Figure  4.9.  In  this  picture  the 
indentation  has  shrunk  down  to  a  minimum  of  about  1". 


Figure  4.7:  Picture  of  indentation  created  by  a  1/16"  nozzle  jet  with  an  exit  velocity 
of  220  ft/s.  The  indentation  is  at  a  relative  maximum  of  about  3/8". 
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Figure  4.8:  Picture  of  indentation  created  by  a  1/16"  nozzle  jet  with  an  exit  velocity 
of  220  ft/s  taken  5  frames  after  Figure  4.7.  The  indentation  is  at  a  relative  minimum 
of  about  1/4". 


Figure  4.9:  Picture  of  indentation  created  by  a  1/8"  nozzle  jet  with  an  exit  velocity 
of  220  ft/s.  The  indentation  is  at  a  relative  maximum  of  about  1.25". 
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Figure  4.10:  Picture  of  indentation  created  by  a  1/8"  nozzle  jet  with  an  exit  velocity 
of  220  ft/s  taken  8  frames  after  Figure  4.9.  The  indentation  is  at  a  relative  minimum 
of  about  1". 

4.3.1  Relationship  between  jet  velocity  and  acoustic  power 

The  data  from  the  experiments  has  been  plotted  to  check  the  dependence  of  the 
radiated  sound  on  the  jet  velocity.  Analytical  predictions  for  a  gas  jet  impinging 
on  a  water  interface  give  a  U 6  dependence  between  jet  exit  velocity  and  acoustic 
pressure  spectrum  in  the  water  [1,2], 

P acoustic  ~  P()U  M  (4.3.1) 

where  M  is  the  jet  Mach  number  relative  to  the  sound  speed  in  water.  If  we 
put  p  =  all13,  then 


20logwp  =  2Dlogwa  +  /320logi0U 


(4.3.2) 
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Plotting  the  pressure  spectrum  as  a  function  of  the  jet  velocity  on  a  logarithmic 
scale  shows  the  slope  j3.  A  slope  of  f3  =  3  would  correspond  to  a  U 6  velocity 
dependence.  Figure  4.11  is  a  plot  of  the  average  acoustic  power  spectrum  vs  the  log 
of  the  velocity  for  the  1/16"  nozzle.  Unfortunately  in  this  case,  the  acoustic  power 
level  does  not  show  any  dependence  on  the  exit  velocity  of  the  gas  jet.  Figure  4.12 
is  a  plot  of  the  average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the 
1/8"  nozzle.  The  slope  of  this  curve  on  a  log  scale  is  i3  —  1.98.  Figure  4.13  is  a 
plot  of  the  average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  3/16" 
nozzle.  The  slope  of  this  curve  on  a  log  scale  is  j3  —  2.55.  Figure  4.14  is  a  plot  of 
the  average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  0.355"  nozzle. 
The  slope  of  this  curve  on  a  log  scale  is  f3  =  4.622.  It  appears  that  there  is  a  relation 
between  the  area  of  the  jet  nozzle  and  the  associated  power  law  for  velocity.  As  the 
diameter  of  the  jet  nozzles  increases,  the  slope  of  the  fitted  curves  /3  increase. 


Average  Pressure  Spectrum  Level  vs  Jet  Exit  Velocity,  2012  D  =  1/16" 
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Figure  4.11:  Average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  1/16" 
nozzle.  This  plot  shows  no  dependence  on  exit  velocity. 
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Average  Pressure  Spectrum  Level  vs  Jet  Exit  Velocity,  2012  D  =  1/8" 
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Figure  4.12:  Average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  1/8" 
nozzle.  The  slope  of  the  data  fit  is  (3  =  1.98. 


Average  Pressure  Spectrum  Level  vs  Jet  Exit  Velocity,  2012  D  =  3/16" 
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Figure  4.13:  Average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  3/16" 
nozzle.  The  slope  of  the  data  fit  is  (3  —  2.55. 
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Average  Pressure  Spectrum  Level  vs  Jet  Exit  Velocity,  2012  D  =  0.355" 


Figure  4.14:  Average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  0.355" 
nozzle.  The  slope  of  the  data  fit  is  3  —  4.622. 


Additional  acoustic  measurements  were  taken  in  August  of  2013.  These 
measurements  were  related  to  a  deflector  plate  experiment  that  ended  up  being 
outside  of  the  scope  of  this  thesis.  But  for  a  couple  of  runs,  the  test  setup  was 
virtually  identical  to  the  measurements  taken  in  2012.  The  3/16"  nozzle  was  the 
only  jet  examined  in  2013.  The  difference  between  these  two  experiments  was  that 
the  nozzle  was  positioned  5.5"  above  the  water.  Figure  4.15  compares  the  data  from 
2012  and  2013  with  an  exit  velocity  of  140  ft/s.  The  spectrums  shows  some  similar 
characteristics  between  the  two  measurements,  especially  the  tone  at  7  kHz.  The 
overall  power  levels  are  considerably  lower  due  to  the  increased  distance  of  the  jet 
from  the  air-water  interface.  Figure  4.16  compares  the  data  from  2012  and  2013 
with  an  exit  velocity  of  220  ft/s.  Again,  the  spectra  show  similar  features  including 
the  tone  at  7  kHz. 
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Hydrophone  Pressure  Power  Spectral  Density,  D  =  3/16"  V  =  140,  L  =  1.75",  5.5" 
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Figure  4.15:  Comparison  of  measured  acoustic  spectra  between  2012  and  2013  with 
jet  exit  velocities  of  140  ft/s. 


Hydrophone  Pressure  Power  Spectral  Density,  D  =  3/16"  V  =  220,  L  =  1.75",  5.5" 


Figure  4.16:  Comparison  of  measured  acoustic  spectra  between  2012  and  2013  with 
jet  exit  velocities  of  220  ft/s. 


It  looks  like  the  power  levels  from  2013  increased  more  rapidly  than  the  power 


levels  in  2012  with  an  increase  in  velocity.  To  check  this,  the  power  law  for  the 
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2013  data  is  plotted  in  Figure  4.17.  The  slope  of  the  curve  fit  for  the  2013  data  was 
/3  =  3.45.  This  is  a  bit  higher  than  the  power  law  predicted  in  2012  for  the  same 
size  nozzle.  The  longer  jet  length  used  in  2013  increases  the  impact  area  of  the  gas 
jet.  This  provides  further  evidence  that  the  power  law  relation  for  velocity  increases 
with  an  increase  in  jet  impact  area. 


Average  Pressure  Spectrum  Level  vs  Jet  Exit  Velocity,  2013 
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Figure  4.17:  Average  acoustic  power  spectrum  vs  the  log  of  the  velocity  for  the  3/16" 
nozzle  taken  in  2013.  The  slope  of  the  data  fit  is  f3  =  3.45. 

The  previous  charts  showed  power  law  calculations  using  the  average  pressure 
level  of  the  entire  frequency  range.  In  order  to  get  an  idea  of  what  was  happening  in 
different  bands  of  frequency,  averages  were  taken  of  the  pressure  levels  for  specific 
ranges  of  frequencies.  The  low,  middle,  and  high  frequency  bins  averaged  the 
spectrum  pressure  level  between  2-7  kHz  7-20kHz,  and  20-50  kHz  respectively. 

Figures  4.18-4.20  show  the  power  law  calculations  for  the  1/8"  nozzle  over  the 
low,  middle,  and  high  frequency  bins.  For  the  1/8"  nozzle  the  low  frequency  bin 
calculated  /3  was  1.04,  the  middle  frequency  calculated  f3  was  1.44,  and  the  high 
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frequency  calculated  f3  was  2.56.  Figures  4.21-4.23  show  the  power  law  calculations 
for  the  3/16"  nozzle  over  the  low,  middle,  and  high  frequency  bins  respectively.  For 
the  3/16"  nozzle  the  low  frequency  bin  calculated  (3  was  2.27,  the  middle  frequency 
bin  calculated  (3  was  2.82,  and  the  high  frequency  bin  calculated  f3  was  2.52.  Figures 
4.24-4.26  show  the  power  law  calculations  for  the  0.355"  nozzle  over  the  low,  middle, 
and  high  frequency  bins  respectively.  For  the  0.355"  nozzle  the  low  frequency  bin 
calculated  j3  was  2.54,  the  middle  frequency  bin  calculated  (3  was  2.75,  and  the  high 
frequency  bin  calculated  f3  was  5.46.  The  power  law  calculations  by  frequency  bin 
are  similar  to  those  calculated  using  the  entire  frequency  range. 


Average  Pressure  Spectrum  in  2-7  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  1/8' 
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Figure  4.18:  Average  acoustic  power  spectrum  between  2-7  kHz  vs  the  log  of  the 
velocity  for  the  1/8"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  (3  =  1.04. 
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Average  Pressure  Spectrum  in  7-20  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  1/8 


Figure  4.19:  Average  acoustic  power  spectrum  between  7-20  kHz  vs  the  log  of  the 
velocity  for  the  1/8"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  /3  =  1.45. 


Average  Pressure  Spectrum  in  20-50  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  1/1 


Figure  4.20:  Average  acoustic  power  spectrum  between  20-50  kHz  vs  the  log  of  the 
velocity  for  the  1/8"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  /3  =  2.56. 
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Average  Pressure  Spectrum  in  2-7  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  3/16 


Figure  4.21:  Average  acoustic  power  spectrum  between  2-7  kHz  vs  the  log  of  the 
velocity  for  the  3/16"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  d  —  2.27. 


Average  Pressure  Spectrum  in  7-20  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  3/1 1 
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Figure  4.22:  Average  acoustic  power  spectrum  between  7-20  kHz  vs  the  log  of  the 
velocity  for  the  3/16"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  j3  =  2.82. 
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Average  Pressure  Spectrum  in  20-50  kHz  Range  Level  vs  Jet  Exit  Velocity,  2012  D  =  3/1 


Figure  4.23:  Average  acoustic  power  spectrum  between  20-50  kHz  vs  the  log  of  the 
velocity  for  the  3/16"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  /3  —  2.52. 


Average  Pressure  Spectrum  Level  in  2-7  kHz  Range  vs  Jet  Exit  Velocity,  2012  D  =  0.35 


Figure  4.24:  Average  acoustic  power  spectrum  between  2-7  kHz  vs  the  log  of  the 
velocity  for  the  0.355"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  [3  —  2.54. 
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Average  Pressure  Spectrum  Level  in  7-20  kHz  Range  vs  Jet  Exit  Velocity,  2012  D  =  0.3E 


Figure  4.25:  Average  acoustic  power  spectrum  between  7-20  kHz  vs  the  log  of  the 
velocity  for  the  0.355"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  (3  =  2.75. 


Average  Pressure  Spectrum  Level  in  20-50  kHz  Range  vs  Jet  Exit  Velocity,  2012  D 


=  0.3 


Figure  4.26:  Average  acoustic  power  spectrum  between  20-50  kHz  vs  the  log  of  the 
velocity  for  the  0.355"  nozzle  taken  in  2012.  The  slope  of  the  data  fit  is  (3  —  5.46. 


68 


4.4  Conclusions 

The  recordings  of  the  underwater  acoustic  pressure  spectrums  for  various  jets 
at  different  velocities  provided  additional  insight  into  the  radiated  sound  from 
supercavitating  vehicles.  The  characteristic  7  kHz  tone  of  the  gas  injection  system 
was  strong  enough  to  radiate  into  the  water.  This  indicates  that  it  is  vital  to 
understand  the  in  air  sound  produced  by  the  ventilation  system  of  a  supercavitating 
vehicle. 

The  relationship  between  jet  exit  velocity  and  radiated  sound  was  also 
explored.  With  the  exception  of  the  1/16"  nozzle,  the  measured  power  law 
relation  was  close  to  the  f/6  relation  predicted  analytically.  The  power  law  relation 
for  the  gas  jets  also  appeared  to  be  related  to  the  jet  impact  area.  The  levels 
of  radiated  sound  from  gas  jets  with  larger  impact  areas  change  more  rapidly 
than  gas  jets  with  smaller  impact  areas.  This  was  particularly  evident  when 
comparing  the  data  from  2012  and  2013.  The  only  difference  between  the  two 
experiments  was  the  distance  the  nozzles  were  positioned  above  the  air-water 
interface.  The  jet  was  positioned  higher  above  the  interface  in  2013,  which  created 
a  larger  jet  impact  area  at  the  surface.  The  power  law  related  to  exit  velocity 
determined  for  the  2013  acoustic  spectrum  levels  was  higher  than  the  data  from  2012. 


Chapter  5 


Acoustic  measurements  of  gas  jets 
generating  bubble  entrainment  at  the 
air-water  interface 

This  chapter  involves  the  acoustic  measurements  of  the  underwater  noise  from 
gas  jets  impinging  on  the  air-water  interface  with  bubble  entrainment.  Section 
5.1  describes  the  experimental  setup  and  procedure  for  taking  the  measurements. 
Section  5.2  contains  descriptions  of  the  data  analysis,  including  any  post-processing, 
and  a  discussion  of  sources  of  uncertainty  in  the  data.  Section  5.3  presents  the 
results  in  graphical  format  and  discusses  some  of  the  trends  in  the  data.  Section  5.4 
compares  the  measured  data  from  the  experiment  to  predicted  sound  levels  from 
the  analytical  model  of  sound  traveling  through  a  bubble  hemisphere.  Section  5.5 
summarizes  the  important  conclusions  drawn  from  the  experiments. 


5.1  Experimental  setup  and  procedures 

The  setup  and  procedure  for  this  measurement  was  very  similar  to  the  methods  used 
to  measure  the  underwater  sound  from  the  gas  jets  without  bubble  entrainment. 
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Experiments  were  performed  at  the  Acoustic  Test  Facility  at  NUWC  to  measure  the 
noise  generated  by  gas  jets  hitting  the  air-water  interface  with  bubble  entrainment 
on  Aug.  20-21  of  2012  at  the  Acoustic  Test  Facility  (ATF)  in  NUWC  Newport,  HI. 
The  ATF  is  a  large  indoor  laboratory  that  specializes  in  accurate  underwater  sound 
recordings.  The  facility  features  a  700,000  gallon  test  tank,  with  dimensions  of  60 
ft.  long,  40  ft.  wide,  and  35  ft.  deep.  Portions  of  the  tank  are  lined  with  acoustic 
wedges  to  dampen  out  reflections. 

The  gas  ventilation  system  was  positioned  on  a  catwalk  above  the  acoustic 
tank.  The  jet  nozzle  was  attached  to  a  lowering  arm  and  positioned  the  jet  1.75" 
above  the  interface.  A  total  of  4  hydrophones  were  positioned  in  the  tank  beneath 
the  catwalk  to  record  the  underwater  sound  as  shown  in  Figure  4.1.  The  primary 
hydrophone  was  oriented  horizontally  and  positioned  3  m  directly  below  the  gas 
jet.  Another  2  hydrophones  were  oriented  vertically  and  positioned  1  m  away  from 
the  jet  and  1  m  below  the  surface.  The  fourth  hydrophone  was  again  oriented 
horizontally  and  positioned  1  m  away  from  the  jet  and  1  m  below  the  surface. 
The  orientation  of  the  hydrophones  is  important  because  of  the  directivity  of  the 
hydrophones.  The  hydrophones  are  least  sensitive  in  the  direction  towards  the  cable 
they  are  attached  to.  The  fourth  hydrophone  was  added  so  there  would  be  two 
phones  oriented  in  the  same  direction  for  comparison. 

In  addition  to  the  hydrophones,  an  underwater  camera  was  used  to  take  images 
of  the  bubble  clouds  generated  by  the  impinging  jets.  In  order  to  get  better  contrast 
in  the  underwater  images,  a  black  background  was  submerged  behind  the  region 
where  the  jets  impacted  the  surface  of  the  water.  Underwater  lights  were  positioned 
off  to  the  side,  projecting  light  towards  the  bubble  cloud  and  the  black  background. 
Finally,  a  meter  stick  was  submerged  near  the  entrained  bubble  cloud  to  provide  a 
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sense  of  scale  to  the  images. 

Data  was  acquired  from  the  hydrophones  with  the  data  acquisition  system  at 
the  ATF.  The  acoustic  data  was  sampled  at  a  rate  of  150  kHz.  The  time  series  was 
then  converted  to  a  power  spectrum  using  16,384  points,  and  the  final  spectrum 
recorded  was  a  result  of  16  spectral  averages.  For  the  bubble  entrainment  cases,  the 
1/16"  and  1/8"  nozzles  were  used.  It  was  observed  in  preliminary  trials  that  the 
3/16"  and  .355"  nozzles  were  not  effective  at  entraining  bubbles.  The  nozzles  with 
larger  diameters  tend  to  produce  relatively  stable  dimples  at  the  surface,  with  few 
bubbles.  It  is  easier  for  the  jets  with  smaller  diameters  to  break  through  the  surface 
tension  to  entrain  bubbles  beneath  the  surface. 

For  each  trial,  the  following  procedure  was  used  to  acquire  the  data.  First, 
recording  of  the  gas  instrumentation  sensors  was  started  in  Labview.  Then  the  gas 
flow  was  turned  on,  and  the  pressure  regulator  was  adjusted  to  get  the  gas  to  the 
desired  flow  rate.  This  flow  rate  was  estimated  using  incompressible  flow  theory, 
along  with  the  nozzle  diameter  and  desired  exit  velocity.  After  the  desired  flow  rate 
was  achieved,  the  acoustic  data  acquisition  system  was  started.  After  the  system 
had  performed  the  16  spectral  averages,  gas  flow  was  shut  off,  followed  by  the  gas 
data  acquisition  system. 


5.2  Data  analysis 

5.2.1  Post-processing 


Data  was  acquired  from  the  hydrophones  with  the  data  acquisition  system  at  the 
ATF.  The  acoustic  data  was  sampled  at  a  rate  of  150  kHz.  The  time  series  was  then 
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converted  to  a  power  spectrum  using  16,384  points,  and  the  final  spectrum  recorded 
was  a  result  of  16  spectral  averages.  The  data  was  then  smoothed  over  frequency 
using  a  running  average  with  a  frequency  window  of  82  Hz  and  90  percent  overlap. 

5.2.2  Uncertainty 

The  principal  sources  of  uncertainty  in  this  experiment  are  from  the  chaotic 
interaction  of  the  high  speed  gas  jets  with  the  air-water  interface.  The  entrained 
bubbles  beneath  the  surface  generally  formed  an  ellipsoid  shape,  but  the  precise  size 
of  the  cloud  was  rapidly  changing.  The  size  and  distribution  of  the  bubbles  withing 
the  cloud  were  also  rapidly  changing.  The  high  impact  jets  also  splashed  quite  a 
bit  of  water  into  the  air.  These  drops  fell  back  down  to  the  interface  similar  to 
raindrops,  and  could  have  impacted  the  acoustic  measurements.  Measurements  of 
the  exit  velocities  for  these  high  speed  jets  indicate  that  the  gas  jet  injection  system 
can  successfully  maintain  a  steady  exit  velocity  for  these  high  flow  rates.  The  only 
other  source  of  error  could  be  from  the  hydrophones.  The  acoustic  data  is  converted 
to  decibels,  and  due  to  their  logarithmic  nature  the  error  is  customarily  marked  as 
+/  -  1  dB. 


5.3  Results 

The  data  from  the  bubble  entrainment  cases  has  been  plotted  in  Figures  5. 1-5.2. 
Figure  5.1  shows  the  smoothed  pressure  data  for  the  1/16"  jet  with  the  pressure 
regulator  set  to  50  psi.  A  plume  of  bubbles  was  created  that  extended  about  4" 
below  the  surface  of  the  water,  as  seen  in  Figure  5.3.  Figure  5.2  shows  the  smoothed 


73 


pressure  data  for  the  1/8"  jet  also  with  the  pressure  regulator  set  to  50  psi.  Figure 
5.4  shows  a  plume  of  bubbles  was  created  that  extended  about  6.5"  below  the 
surface  for  the  1/8"  nozzle.  The  recorded  sound  spectrum  of  the  1/16"  jet  and  the 
1/8"  jet  are  of  similar  magnitude  despite  the  larger  impact  area  created  by  the  1/8" 
jet. 


Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/16” 
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Figure  5.1:  Pressure  spectrum  for  1/16"  nozzle,  pressure  regulator  set  to  50  psi  for 
bubble  entrainment. 
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Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/8" 


Figure  5.2:  Pressure  spectrum  for  1/8"  nozzle,  pressure  regulator  set  to  50  psi  for 
bubble  entrainment. 


Figure  5.3:  Underwater  image  of  the  bubble  cloud  produced  by  the  1/16"  nozzle  set 
to  50  psi  for  bubble  entrainment. 
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Figure  5.4:  Underwater  image  of  the  bubble  cloud  produced  by  the  1/8"  nozzle  set 
to  50  psi  for  bubble  entrainment. 

5.4  Conclusions 

When  the  pressure  is  set  to  50  psi,  the  gas  jet  injection  system  is  capable  of 
entraining  bubbles  beneath  the  air- water  interface  with  the  1/8"  and  1/16" 
nozzles.  The  measurements  of  sound  traveling  through  a  bubble  hemisphere 
provided  experimental  evidence  of  an  acoustic  signal  being  amplified  over  a  range 
of  frequencies  by  the  bubble  hemisphere.  Observed  dimensions  of  the  bubble 
hemispheres  will  be  used  as  an  input  to  the  theoretical  model  developed  later  in 


6.3.5. 


Chapter  6 


Sound  generated  by  a  gas-jet 
impingement  on  a  bubbly  gas-water 
interface 

This  chapter  describes  the  predicted  sound  radiation  when  a  gas  jet  impacts  a  bub¬ 
bly,  two-phase  interface.  Section  6.1  presents  the  physical  model  and  assumptions 
made  in  the  initial  formulation  of  the  model.  In  the  absence  of  bubbles,  Foley  et 
al  [1]  showed  that  the  radiated  sound  was  equivalent  to  that  produced  by  a  dipole 
source  whose  axis  is  normal  to  the  interface  and  that  the  amplitude  is  dependent  on 
the  unsteady  force  applied  by  the  gas  jet.  Section  6.2  outlines  the  important  steps 
behind  Foley’s  derivation.  Section  6.3  describes  the  influence  on  Foley’s  result  when 
a  cloud  of  bubbles  is  created  beneath  the  jet  impact  point.  Section  6.4  discusses  the 
effect  of  a  bubbly  layer  at  the  air-water  interface  on  the  radiated  sound.  Section 
6.5  summarizes  the  previous  sections  and  discusses  the  implications  of  the  predicted 
results. 
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6.1  The  acoustic  problem 

Foley  et  al.  [1,2]  assumed  that  the  unsteady  forces  produced  by  the  20  jets  on  the 
gas-water  interface  of  the  experimental  supercavity  can  be  regarded  as  statistically 
independent.  Each  jet  is  turbulent  and  impinges  at  near  normal  incidence,  and 
together  they  form  a  set  of  uncorrelated  acoustic  sources.  The  spectrum  of  the 
overall  far  held  jet  noise  can  therefore  be  obtained  by  simple  addition  of  the  noise 
spectra  of  the  individual  jets. 


Figure  6.1:  Foley’s  model  used  to  estimate  the  sound  generated  by  an  unsteady  jet 
incident  normally  on  the  gas-water  interface. 

6.2  Foley’s  problem 

In  practice  the  diameter  D  of  the  jet  impact  area  on  the  interface  is  much  smaller 
than  the  local  cavity  radius,  so  that  the  interface  can  be  assumed  locally  plane. 
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Foley  et  al.  [1,2]  confined  attention  to  those  components  of  the  sound  having 
wavelengths  (in  the  water)  larger  than  D.  This  permits  the  impact  region  to  be 
reduced  to  a  fluctuating  force  applied  at  a  single  point  on  the  interface.  The  model 
of  sound  generation  at  high  frequencies  is  illustrated  schematically  in  Figure  6.1. 
The  locally  plane  jet-impact  region  of  the  gas-water  interface  is  taken  to  coincide 
with  the  plane  £3  =  0  of  the  rectangular  coordinate  system  x  =  (xi,X2,x3),  where 
the  water  of  density  pQ  and  sound  speed  ca  lies  in  x3  >  0.  The  Mach  number  of  the 
water  flow  relative  to  the  cavity  is  assumed  to  be  small  enough  that  convection  and 
refraction  of  sound  by  the  mean  flow  can  be  neglected,  such  that  small  amplitude 
acoustic  pressures  p(x,  t)  in  the  water  satisfy 

(^-V2)p  =  0.  (6.2.1) 

The  idealized  gas  jet  has  circular  cross-section,  is  coaxial  with  the  x3  axis 
and  impinges  normally  on  the  interface  from  within  the  cavity  (x:i  <  0).  The  jet 
exhausts  at  mean  speed  Ua  from  a  circular  nozzle  of  diameter  Da  situated  a  distance 
t  from  the  interface  (where  D0  <C  £),  and  is  assumed  to  exert  an  unsteady  force 
on  the  water  over  a  circular  region  of  diameter  D  <C  t.  If  the  dominant  acoustic 
wavelengths  produced  by  the  jet  are  much  larger  than  D  the  distributed  force  may 
be  regarded  as  concentrated  at  a  point  on  the  interface,  corresponding  to  a  singular 
distribution  of  surface  pressure  p  =  FS(xi)S(x2),  where  F  is  the  net  force.  The  force 
consists  of  a  mean  component  together  with  a  fluctuating  part;  only  the  unsteady 
part  of  the  force  is  responsible  for  sound  production,  and  we  shall  henceforth  denote 
it  by  F(t)  and  assume  it  to  be  a  stationary  random  function  of  the  time.  The 
unsteady  pressure  elsewhere  on  the  interface  can  be  neglected  because  of  the  very 
large  differences  in  the  mean  gas  and  water  densities. 
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The  mathematical  problem  of  sound  production  by  a  single  jet  is  accordingly 
governed  by  equation  (6.2.1)  in  x3  >  0  subject  to  the  condition  p  =  F(t)5(xi)5(x2) 
on  £3  =  0.  The  solution  with  outgoing  wave  behavior  is  the  dipole  pressure  field  [8, 
9] 


P(x,t)  = 


r-v_/ 


d  I"  1  /  Ixl 

—  -  - Fit  —  — 

dx3  L2vr|x|  V  c0 

cos  6  d F  (  |x|\ 

27tc0|x|  dt  \  cQ  J  ’ 


x 


,  X3  >  0, 

— »  00,  cos  6 


Ixl 


(6.2.2) 

(6.2.3) 


The  frequency  spectrum  <h^(a;)  of  the  jet  force  is  defined  such  that 


(■ F2(t ))  =  /  $FF(cu)dcu, 
Jo 


(6.2.4) 


where  the  angle  brackets  (  )  denote  an  ensemble  or  time  average.  The  corresponding 
acoustic  pressure  frequency  spectrum  <f>pp(x,  uj)  satisfies 


(p2(x,t)}  =  /  $pp(x,uj)duj 

Jo 


(6.2.5) 


and  it  follows  from  (6.2.3)  that 


cos^  0 

4>pp(x,cu)  «  27r2|x|2  ko®FF(u)’  lxl^°°- 


(6.2.6) 


where  kQ  =  oj/c0  is  the  acoustic  wavenumber. 


6.2.1  Comparison  of  acoustic  model  to  measured  results 


Foley’s  model  was  the  most  appropriate  to  use  for  cases  where  the  interface  was 
relatively  smooth  with  no  bubbles.  The  measured  unsteady  force  spectra  described 
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in  Chapter  3  were  used  as  inputs  to  Foley’s  original  model.  The  measurements  from 
each  force  sensor  were  added  together  for  a  cumulative  force  input  for  the  model  as 
shown  in  Equation  6.2.7. 

FtOTAL  =  Fai\  +  FC528  +  Fq  501  +  Fq 502  +  Fq 532  +  Fq  530  +  ^519  (6.2.7) 

The  total  force  was  calculated  by  summing  up  the  contributions  from  the  individual 
force  sensors  in  Newtons.  The  predicted  acoustic  model  for  the  1/16"  nozzle  is 
shown  in  Figure  6.2.  For  this  nozzle,  the  force  data  hits  the  noise  floor  at  about  8 
kHz.  The  predicted  models  for  the  1/8",  3/16",  and  0.355"  are  shown  in  Figures 
6. 3-6. 5.  The  spikes  in  the  model  at  32  kHz  are  from  the  resonant  frequency  of  the 
force  sensors. 

Figures  6. 6-6. 8  shows  a  comparison  of  the  measured  and  modeled  acoustic 
spectrums  for  the  1/8",  3/16",  and  0.355"  nozzles  respectively  with  an  exit  velocity 
of  220  ft/s.  The  predicted  levels  in  all  of  the  models  were  about  15  dB  below  the 
measured  acoustic  power  level.  There  are  a  few  possible  explanations  for  the  under 
prediction  of  the  current  acoustic  model.  It  was  assumed  during  the  derivation  of 
the  model  that  the  area  of  jet  impact  could  be  represented  with  a  delta  function.  In 
reality,  the  jet  exerts  force  on  a  small  circular  area  on  the  interface.  There  could  be 
additional  sound  radiated  into  the  water  that  is  not  accounted  for  due  to  the  use 
of  a  delta  function  to  represent  the  impact  area.  Also,  the  current  model  does  not 
account  for  any  in-air  jet  noise  transmitting  through  the  interface.  Sound  of  the  air 
jet  through  the  interface  is  not  the  primary  mechanism  of  radiated  sound  in  the 
water,  but  it  could  contribute  some  to  the  overall  energy. 
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Acoustic  Model  of  Pressure  Spectrum  for  1/16"  Nozzle 


Figure  6.2:  Predicted  acoustic  model  for  the  1/16"  nozzle. 


Acoustic  Model  of  Pressure  Spectrum  for  1/8"  Nozzle 


Figure  6.3:  Predicted  acoustic  model  for  the  1/8"  nozzle. 
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Acoustic  Model  of  Pressure  Spectrum  for  3/16"  Nozzle 


Figure  6.4:  Predicted  acoustic  model  for  the  3/16"  nozzle. 


Acoustic  Model  of  Pressure  Spectrum  for  0.355"  Nozzle 


Figure  6.5:  Predicted  acoustic  model  for  the  0.355"  nozzle. 
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Measured  and  Modeled  Acoustic  Power  Spectral  Density,  D  =  1/8",  V  =  220  ft/s 


Figure  6.6:  Comparison  of  measured  and  predicted  acoustic  pressure  spectrum  for 
the  1/8"  nozzle  with  an  exit  velocity  of  220  ft/s. 


Measured  and  Modeled  Acoustic  Power  Spectral  Density,  D  =  3/16",  V  =  220  ft/s 


Figure  6.7:  Comparison  of  measured  and  predicted  acoustic  pressure  spectrum  for 
the  3/16"  nozzle  with  an  exit  velocity  of  220  ft/s. 


84 


Measured  and  Modeled  Acoustic  Power  Spectral  Density,  0.355",  V  =  220  ft/s 


Figure  6.8:  Comparison  of  measured  and  predicted  acoustic  pressure  spectrum  for 
the  0.355"  nozzle  with  an  exit  velocity  of  220  ft/s. 

The  relationship  between  jet  exit  velocity  and  the  modeled  acoustic  pressure 
spectrum  level  was  determined  with  the  same  logic  applied  in  section  4.3.1.  The 
pressure  spectrum  level  was  averaged  between  a  range  of  frequencies  and  plotted 
against  the  log  of  the  exit  velocity.  A  curve  fit  yielded  a  coefficient  f3  which 
determined  the  power  law  relation  between  jet  exit  velocity  and  pressure  spectrum 
level.  The  noise  floor  associated  with  the  force  spectrum  measurements  was  used  as 
the  limit  for  the  frequency  range  of  the  analyzed  model  pressure  levels. 

Figure  6.9  shows  the  power  law  relation  between  exit  velocity  and  modeled 
pressure  spectrum  levels  for  the  1/8"  nozzle  using  a  frequency  range  of  2-7  kHz. 
The  model  hits  the  noise  floor  of  the  force  spectrum  measurements  above  7  kHz.  A 
curve  fit  for  the  average  pressure  spectrums  within  the  2-7  kHz  frequency  ranges 
yielded  a  (3  of  1.82. 

For  the  3/16"  nozzles,  the  model  hits  the  noise  floor  at  17  kHz.  The  power 
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law  relation  f3  was  calculated  over  two  frequency  bins,  the  first  ranges  from  2-7  kHz 
while  the  second  ranges  from  7-17  kHz.  Figures  6.10-6.11  plot  the  calculated  curve 
fits  for  the  3/16"  nozzles  over  the  first  and  second  frequency  bins  respectively.  This 
yielded  a  /3  of  2.21  between  2-7  kHz,  and  a  (5  of  2.43  between  7-17  kHz. 

For  the  0.355"  nozzle,  there  were  no  issues  with  the  noise  floor  of  the  force 
spectrum  measurements.  But  the  modeled  data  in  the  upper  frequencies  above 
20  kHz  was  significantly  affected  by  the  internal  resonance  of  the  force  sensors. 
Therefore  modeled  data  above  20  kHz  was  not  used  in  the  power  law  calculations. 
The  power  law  relation  f3  was  calculated  over  two  frequency  bins,  the  first  ranges 
from  2-7  kHz  while  the  second  ranges  from  7-20  kHz.  Figures  6.12-6.13  plot  the 
calculated  curve  fits  for  the  0.3554"  nozzles  over  the  first  and  second  frequency 
bins  respectively.  A  /3  of  2.19  was  calculated  between  2-7  kHz,  and  a  f3  of  3.63  was 
calculated  between  7-20  kHz.  The  power  law  relations  determined  for  the  modeled 
pressure  spectrum  and  jet  exit  velocity  is  similar  in  scale  to  the  /Ts  calculated  for 
the  measured  acoustic  spectrum. 


Modeled  Average  Pressure  Spectrum  in  2-7  kHz  Range  vs  Jet  Exit  Velocity,  D  =  1/8" 


Figure  6.9:  Average  modeled  acoustic  power  spectrum  vs  the  log  of  the  velocity  for 
the  1/8"  nozzle  over  the  2-7  kHz  frequency  range.  The  slope  of  the  data  fit  is  [3  — 


1.82. 


Modeled  Average  Pressure  Spectrum  in  2-7  kHz  Range  vs  Jet  Exit  Velocity,  D  =  3/16" 
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Figure  6.10:  Average  modeled  acoustic  power  spectrum  vs  the  log  of  the  velocity  for 
the  3/16"  nozzle  over  the  2-7  kHz  frequency  range.  The  slope  of  the  data  fit  is  [3  — 
2.21. 
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Modeled  Average  Pressure  Spectrum  in  7-17  kHz  Range  vs  Jet  Exit  Velocity,  D  =  3/16' 


Figure  6.11:  Average  modeled  acoustic  power  spectrum  vs  the  log  of  the  velocity  for 
the  3/16"  nozzle  over  the  7-17  kHz  frequency  range.  The  slope  of  the  data  fit  is  3  — 
2.43. 


Modeled  Average  Pressure  Spectrum  in  2-7  kHz  Range  vs  Jet  Exit  Velocity,  D  =  .355" 
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Figure  6.12:  Average  modeled  acoustic  power  spectrum  vs  the  log  of  the  velocity  for 
the  0.355"  nozzle  over  the  2-7  kHz  frequency  range.  The  slope  of  the  data  fit  is  [3  — 
2.19. 


Modeled  Average  Pressure  Spectrum  in  7-20  kHz  Range  vs  Jet  Exit  Velocity,  D  =  .355' 
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Figure  6.13:  Average  modeled  acoustic  power  spectrum  vs  the  log  of  the  velocity  for 
the  0.355"  nozzle  over  the  7-20  kHz  frequency  range.  The  slope  of  the  data  fit  is  f3 
=  3.63. 


6.3  Influence  of  a  bubbly  cloud 

This  derivation  was  motivated  by  observations  of  gas  jets  impinging  on  the  air-water 
interface  at  high  velocities.  At  sufficient  speeds,  the  gas  jets  would  entrain  a  cloud 
of  bubbles  beneath  the  surface  of  the  water.  These  clouds  were  generally  ellipsoid  in 
shape,  but  the  precise  geometry  of  the  cloud  would  fluctuate  rapidly.  The  best  start 
at  an  analytical  model  for  this  observed  phenomenon  was  to  begin  with  a  bubble 
hemisphere. 


6.3.1  Mean  properties  of  the  cloud 

The  influence  on  sound  production  by  a  jet  of  a  bubbly  interface  between  the 
cavity  and  the  water  will  be  investigated  by  modifying  Foley’s  [1,2]  planar  interface 
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approximation  in  the  manner  illustrated  in  Figure  6.14.  It  is  assumed  that  the 
bubble  cloud  is  a  homogeneous  hemisphere  with  radius  R  from  the  coordinate  origin 
.  It  will  be  assumed  that  compression  waves  within  the  layer  have  wavelengths  that 
are  sufficiently  large  compared  to  the  characteristic  bubble  diameter  that  the  layer 
can  be  modeled  analytically  as  a  homogeneous  acoustic  medium  of  whose  density  p 
and  sound  speed  c  are  given  by  [12,38] 


P  =  (1  -  v)Po  +  °Pg 


c2  ~ 


7 gPo 


pQcr{  1  -  cr) 


(6.3.1) 


*3 

A 

sound 


Figure  6.14:  Production  of  sound  by  a  jet  incident  on  the  water  interface  in  the 
presence  of  a  hemispherical  bubble  cloud. 


where  and  pg  are  the  ratio  of  specific  heats  and  the  mean  density  of  the  gas 
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phase,  which  is  assumed  to  occupy  a  fractional  volume  a  within  the  layer,  and  pa  is 
the  mean  cavity  pressure.  The  precise  value  of  a  has  not  been  measured,  but  it  is 
unlikely  to  exceed  0.1  and  probably  exceeds  0.01,  so  that  typically  c/ca  ~  1/15  1. 

Similarly,  the  specific  acoustic  impedance  Z  of  the  bubbly  layer  (relative  to  water) 

is 

Z  =  — .  (6.3.2) 

PoCo 


where  typically  Z  «  0.07. 


At  high  frequencies  the  effective  wave  speed  c  in  the  cloud  is  dependent  on  the 
value  of  the  mean  resonant  frequencies  of  the  individual  bubbles.  This  is  related  to 
the  mean  bubble  radius  Rb  by  the  approximate  formula  [6,12,30,38-40] 


ub 


3jgPo\  1/2 

poRl ) 


(6.3.3) 


For  acoustic  waves  with  time  dependence  of  the  type  e  lut  we  have  [6,40] 


4~4  +  4[l-^ 


u 


-i 


if) —  ,  a  <C  1. 

Ub, 


(6.3.4) 


The  empirical  damping  coefficient  r\  has  been  estimated  by  Fox  et  al,  [32,40]  to  be 
about  /.  This  value  of  r/  is  dependent  on  the  resonant  frequency  of  the  bubbles. 
The  effects  of  adjusting  the  value  of  r/  will  be  investigated  later  in  6.3.5. 


6.3.2  Calculation  of  the  radiated  sound 


Let 


F{t) 


F(uj)e-iwtdt: 


(6.3.5) 
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where  F(uj)  is  the  Fourier  transform  of  the  force  F(t).  The  pressure  distribution  on 
the  interface  x3  =  0  can  then  be  expressed  by 


p(xi,x2,0,t)  =  F(t)S(x1)S(x2) 


F{uj)5(x1)5(x2)e-luJtduj. 


(6.3.6) 


The  first  step  in  calculating  the  radiated  sound  is  to  solve  for  the  time-harmonic 
pressure  with  amplitude  of  F(u>)5(xi)5(x2)  applied  at  x 3  =  0.  The  radiated  sound 
is  then  obtained  by  multiplying  by  e~lu,t  and  integrating  over  —  00  <  cu  <  00. 

The  main  theme  in  this  calculation  is  considering  boundary  conditions,  and 
ensuring  that  there  is  continuity  of  pressure  and  normal  velocity.  Due  to  symmetry, 
the  pressure  radiated  from  the  interface  in  x3  >  0  depends  only  on  the  radial 
distance  r  =  |x|  and  the  polar  angle  9  measured  from  the  positive  x3  direction.  The 
time-harmonic  solution  for  this  type  of  sound  radiation  is 


P(r,  0,  w)  = 


ik%hi(kar)  cos  6 
2n 


FH, 


0<e<l’ 


(6.3.7) 


where 

d  ( Fx 

^  =  ^(t 

is  the  first  order  spherical  Hankel  function  [12,41],  The  boundary  condition  at 
X3  =  0  is  satisfied  because  cos 9  =  x$/r,  and  as  r  =  {x\  +  x%  +  x\)^  — y  0  the  right 
hand  side  of  equation  (6.3.7)  becomes 


x3F(u) 

2n(xf  +  x%  +  x\ )i 


->•  F(u>)5(xi)5(x2) 


The  next  region  to  consider  is  inside  the  bubble  cloud.  The  time-harmonic 
solution  which  also  satisfies  the  same  boundary  condition  at  £3  =  0  the  same 
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boundary  condition  is 

£^2  £Qg  0  /  ^  _  _  \ 

p(r,6,  uj)  =  — yFhi(kr)  +  Ajx(kr)j  ,  r  <  R,  (6.3.8) 

where  k  =  uj/c  is  the  acoustic  wavenumber  within  the  cloud,  j\{x)  is  the  first  order 
spherical  Bessel  function  [41],  and  the  coefficient  A  =  A{u)  still  needs  to  be  solved 
for. 

In  the  clear  water  beyond  the  bubble,  we  use  a  similar  logic  to  write  the 
time-harmonic  solution 

?  ms  f) 

p(r,0,u>)  —  — - - Bhi(k0r),  r  >  R,  (6.3.9) 

2  7T 

where  the  coefficient  B  =  B(lo )  remains  undetermined. 

The  coefficients  A  and  B  are  solved  from  satisfying  the  continuity  of  pressure 
and  normal  velocity  at  r  =  R.  The  normal  velocity  is  related  to  the  pressure  by 
(l/p)dp/dr,  and  so 


Fh^kR)  +  Aji(kR) 
Fh[(kR)  +  Aj[(kR) 


=  Bhiik0R ), 

=  (Z/f)Bh[(k0R) 


(6.3.10) 


where  the  prime  indicates  differentiation,  and 
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u  .  u 
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(6.3.11) 
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Solving  the  system  of  equations  in  (6.3.10)  for  B  and  using  the  Wronksian 
hi(x)j[(x)  -  h[(x)ji(x)  =  -i/x2  [41],  gives 


B  M 


-iF{u) 

(. kR )2  [h^koRMihR)  -  {Z/f)h'i{k0R)h{kR)\ ' 


(6.3.12) 


This  coefficient  can  then  be  substituted  into  equation  (6.3.9).  The  pressure  outside 
the  cloud  is  then  found  by  multiplying  this  result  by  e~lut  and  integrating  over  all 
frequencies 


cosd  r°°  F((jj)hi(k0\x\)e  lu}tdu: 

P(x’ t]  =  FnFP  MWp)  -  (Z/f^ikoRMkRY 


x|  >  R  (6.3.13) 


6.3.3  The  far  field  pressure 

At  ranges  |x|  that  are  much  greater  than  the  characteristic  acoustic  wavelength,  the 
asymptotic  solution  h\(x)  —  elx /x,  |x|  1,  can  be  used  in  (6.3.13)  which  gives 

,  x  —  cos  8 

79(X  t)  ~  -  /  - — - — -  X  — y  oo 

V  ’  27T^|x|  7.^  ’  1  1 

(6.3.14) 

In  the  absence  of  the  bubble  cloud,  k  — »  k0,  and  the  denominator  of  the 
integrand  in  (6.3.13)  becomes  —i/(k0R)2,  and  formula  (6.3.13)  becomes  the  exact 
result  (6.2.2)  for  clear  water.  The  approximation  in  (6.3.14)  also  reduces  to  the 
clear  water  far  field  approximation  (6.2.3).  The  radiated  acoustic  field  exhibits  a 
dipole  directivity  with  and  without  the  presence  of  the  bubble  cloud. 
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6.3.4  The  acoustic  pressure  frequency  spectrum 


The  definitions  of  the  force  and  pressure  spectrums  in  (6.2.4)  and  (6.2.5)  allows  the 
acoustic  pressure  spectrum  to  be  expressed  as 


<f>pP(x,  u)  « 


cos2  9  k ;l$FF(u) 


47r2|x|2(A;0i?)4  [h^KR^ikR)  -  {Z  /  f^Mj^kR^ 


2  ’ 


x  — >  oo. 


(6.3.15) 


When  the  void  fraction  a  of  the  bubble  cloud  goes  to  zero,  (6.3.15)  becomes 
equivalent  to  the  pressure  spectrum  predicted  by  Foley  [1,2]. 

The  acoustic  power  gain  G(u)  produced  by  the  bubble  cloud  is  determined  by 
the  terms  in  the  denominator  related  to  the  composition  of  the  bubble  cloud 


G(u)  =  —10  log10  {{k0Rf  [h^koR^kR)  -  [Zj f)h[ (k0R)ji(kR)] 2 j  dB.  (6.3.16) 


There  are  peaks  in  the  gain  at  the  resonant  frequencies  of  the  bubble  cloud, 
roughly  corresponding  to  the  interior  oscillations  of  a  hemisphere  of  fluid  with  a  flat 
circular  boundary.  This  can  be  seen  in  Figure  6.15  for  the  case  in  which  the  mean 
bubble  resonant  frequency  fb  =  Ub/Zn  =  50  kHz.  Resonance  frequencies  over  for 
other  values  of  fb  show  similar  characteristics. 


A  more  realistic  model  of  the  acoustic  gain  will  take  into  account  the  time 


dependent 
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Figure  6.15:  Mean  gain  G(u)  in  the  sound  pressure  level  produced  by  a  hemispherical 
bubble  cloud  of  nominal  radius  5  cm  for  different  values  of  the  bubble  resonance 
frequency  fb .  The  ’Raw  gain’  is  calculated  without  averaging  with  respect  to  the 
cloud  radius  when  fb  =  50  kHz. 


state  of  the  composition  of  the  bubble  cloud.  This  is  accomplished  by  assuming  a 
uniform  variation  of  the  bubble  cloud  radius  over  the  range  R.mvn  <  R  <  Rmax ■  The 
mean  gain  can  then  be  expressed  as 


(^H)  =  —10  log10(A) 


(6.3.17) 


where 

_  1  fRrnax  dR 

Rmax  Rmin  jRmiu  (. KR )4  [h^koR^kR)  -  (Z/ f)h'1(k0R)ji(kR)]2 

(6.3.18) 

Figure  6.15  also  shows  the  variation  of  the  mean  gain  (G(uu))  for  a  bubble  cloud 
of  nominal  radius  R  =  5  cm  and  for  fb  =  10,30,50  kHz.  The  radius  of  the  cloud 
fluctuates  between  Rm.in  =  Rj 2  and  Rmax  —  3/2/2.  The  averaging  removes  the  large 
peaks  in  the  gain,  and  tends  to  yield  a  more  uniform  gain  of  20-30  dB  over  a  range 
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of  frequencies  below  the  resonant  frequencies  of  the  bubbles.  At  higher  frequencies, 
the  predicted  gain  is  negative,  and  decreases  rapidly  with  increasing  frequency. 


6.3.5  Comparisons  with  measured  acoustic  power 

The  model  for  acoustic  transmission  through  a  bubble  hemisphere  was  used  to 
predict  any  acoustic  signal  gains  through  the  bubble  hemisphere.  The  hemisphere 
radius  and  bubble  size  inputs  used  in  the  model  were  based  on  the  images  taken 
of  the  entrained  bubbles.  The  bubble  clouds  had  an  ellipsoid  shape,  but  the 
simpler  hemispherical  model  still  provided  some  good  results.  Figure  6.16  shows  the 
predicted  acoustic  gain  through  a  bubble  hemisphere  with  a  radius  of  2.5"  composed 
of  tiny  bubbles  with  individual  radii  of  0.032".  The  thick  blue  curve  shows  the 
result  of  averaging  the  spectrum  for  different  bubble  hemisphere  radii.  The  nominal 
radius  was  2.5"  ±  0.62".  This  helps  to  model  a  hemisphere  which  is  rapidly  changing 
size.  The  model  predicts  an  increase  of  about  15  dB  from  about  800-2000  Hz. 
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Mean  Gain,  Hemisphere  Radius:  2.5",  Bubble  Radius:  0.032" 


Figure  6.16:  Predicted  acoustic  gains  through  a  bubble  hemisphere  with  a  radius  of 
2.5"  composed  of  bubbles  with  individual  radii  of  0.032". 

Figures  6.17-6.18  zooms  in  on  the  0-5  kHz  range  of  the  pressure  spectrums  for 
the  1/16"  and  1/8"  pressure  spectrums.  This  view  makes  it  easier  to  see  if  there 
are  any  bumps  in  the  spectrum  that  match  the  prediction  from  the  model.  Both  of 
these  spectrums  show  an  increase  in  signal  through  the  frequency  range  predicted 
by  the  acoustic  model. 
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Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/16" 


Figure  6.17:  Measured  acoustic  pressure  spectrum,  0-5  kHz,  for  the  1/16"  nozzle  set 
to  50  psi. 

Hydrophone  Pressure  Power  Spectral  Density,  D  =  1/8" 


Figure  6.18:  Measured  acoustic  pressure  spectrum,  0-5  kHz,  for  the  1/8"  nozzle  set 
to  50  psi. 


The  measured  force  spectrums  were  used  as  an  input  to  the  bubble  hemisphere 
model.  The  measurements  from  the  individual  force  sensors  were  added  together  to 


99 


produce  a  cumulative  force  that  was  used  as  an  input  to  the  model  (Equation  6.2.7). 
The  measured  force  spectrums  for  the  50  psi  cases  were  significantly  affected  by  the 
resonance  frequency  of  the  force  sensors.  For  this  reason,  the  force  spectrum  was 
corrected  to  match  the  trend  of  the  force  spectrum  data  before  the  resonant  spike 
of  the  force  sensors.  This  straight  line  correction  explains  the  smoothness  of  the 
modeled  curves  at  the  upper  frequencies.  Figure  6.19  plots  the  comparison  of  the 
bubble  hemisphere  acoustic  model  with  the  measured  acoustic  pressure  spectrum 
for  the  1/16"  nozzle  set  at  50  psi.  The  hemisphere  was  modeled  as  having  a  radius 
of  1.18"  composed  of  bubbles  with  internal  radii  of  .0354".  The  effects  of  varying 
the  damping  coefficient  rj  is  explored  in  Figures  6.19-6.20.  A  damping  coefficient  rj 
of  0.1  yields  the  highest  spike  from  the  bubble  hemisphere,  but  also  predicts  a  large 
unrealistic  dip  in  the  predicted  curve.  Larger  values  of  rj  smooth  this  spike  out. 
Figure  6.20  plots  the  comparison  of  the  bubble  hemisphere  acoustic  model  with  the 
measured  acoustic  pressure  spectrum  for  the  1/8"  nozzle  set  at  50  psi.  The  same 
bubble  hemisphere  parameters  were  used  in  this  figure. 
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s  Hemisphere  Model  Compared  to  Measured  Acoustic  Spectrum,  D  =  1/16" 


</) 

<D  50 


-  Measured  Acoustic  Spectrum 

-  Hemisphere  Model,  i]  =  0.5 

-  Hemisphere  Model,  ti  =  0.1 
Hemisphere  Model,  ti  =2 


Figure  6.19:  Comparison  of  modeled  bubble  hemisphere  acoustic  results  and  mea¬ 
sured  pressure  spectrums  for  the  1/16"  set  at  50  psi.  The  hemisphere  was  modeled 
with  a  radius  of  1.18"  composed  of  bubbles  with  internal  radii  of  .0354" 


Bubble  Hemisphere  Model  Compared  to  Measured  Acoustic  Spectrum,  D  =  1/8" 


Figure  6.20:  Comparison  of  modeled  bubble  hemisphere  acoustic  results  and  mea¬ 
sured  pressure  spectrums  for  the  1/8"  set  at  50  psi.  The  hemisphere  was  modeled 
with  a  radius  of  1.18"  composed  of  bubbles  with  internal  radii  of  .0354" 
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6.4  Influence  of  a  bubbly  layer 

This  derivation  was  motivated  by  observations  of  a  supercavitating  vehicle  in  a 
water  tunnel  with  a  rough  and  bubbly  interface  as  seen  in  Figure  1.2.  This  was  a 
result  of  the  ventilating  jets  being  oriented  normal  to  the  interface  instead  of  being 
deflected  to  hit  the  interface  at  a  shallow  angle.  In  a  dynamic  and  unpredictable 
environment  such  as  the  ocean,  it  is  likely  that  the  interface  of  a  supercavity  will 
be  as  perfectly  smooth  as  one  generated  in  a  water  tunnel.  Therefore  the  potential 
presence  and  impact  of  a  bubbly  layer  near  the  interface  is  an  important  case  to 
consider  in  the  application  of  supercavitating  vehicles.  Unfortunately  it  was  not 
within  the  scope  of  this  study  to  experimentally  replicate  a  bubble  layer  near  an 
air-water  interface.  All  of  our  studies  involved  gas  jets  hitting  a  still  body  of  water. 
At  low  impact  velocities  the  gas  jet  would  create  an  indentation  on  the  interface. 
If  the  water  were  moving,  some  of  these  indentations  would  shear  off  and  create  a 
layer  of  bubbles  near  the  interface. 


6.4.1  Mean  properties  of  the  bubble  layer 

The  influence  on  sound  production  by  a  jet  of  a  bubbly  interface  between  the 
cavity  and  the  water  will  be  investigated  by  modifying  Foley’s  [1,2]  planar  interface 
approximation  in  the  manner  illustrated  in  Figure  6.21.  The  bubbles  form  a 
homogeneous  layer  of  nominal  thickness  h  between  the  gas  cavity  and  the  water 
(within  the  interval  0  <  <  h).  It  will  be  assumed  that  compression  waves  within 

the  layer  have  wavelengths  that  are  sufficiently  large  compared  to  the  characteristic 
bubble  diameter  that  the  layer  can  be  modeled  analytically  as  a  homogeneous 
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acoustic  medium  of  whose  mean  density  p  and  sound  speed  c  are  given  by  [12,38] 
Equation  6.3.1. 


Figure  6.21:  Model  used  to  calculate  the  sound  generated  when  the  jet  is  incident 
normally  on  the  lower  face  x%  =  0  of  a  uniform  bubbly  interface  between  the  cavity 
in  x%  <  0  and  water  in  a;3  >  h. 

The  actual  wave  speed  c  in  the  bubble  layer  is  affected  by  the  mean  resonance 
frequency  of  the  bubbles,  Ub,  which  in  turn  depends  on  the  bubble  radius  Rb,  which 
is  given  approximately  by  [6,12,38,40,42]  Equation  6.3.3.  In  a  first  approximation 
when  uCl,  and  for  waves  of  time  dependence  oc  e~lu}t,  we  have  [6,40]  Equation 
6.3.4.  In  terms  of  these  variables,  it  is  convenient  to  introduce  the  specific  acoustic 
impedance  Z  of  the  bubbly  layer  (relative  to  water)  is  given  in  Equation  6.3.2. 
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Then 

F{t)S(Xl)S(x2)  =  JJ°°  F^e^-^^kdu,  (6.4.2) 

where  k  =  (£7,  k2l  0)  and  d2k  =  dkidk2.  The  sound  produced  by  jet  impact  can  now 
be  determined  by  first  considering  that  radiated  into  the  water  when  the  elementary 
pressure  distribution  F(o;)e*l'k'x_aj*')  is  applied  at  x 3  =  0,  over  the  lower  face  of 
the  bubble  layer.  The  aggregate  sound  generated  by  the  jet  is  then  obtained  by 
application  of  the  integral  operator 

dd3  //j')dW  (6'43) 

To  do  this,  the  pressure  field  above  the  interface  (in  £3  >  h)  in  the  simplified 
problem  is  put  in  the  form  [43] 

p  =  Aei{  k-x+7(*3-/0-^  (6.4.4) 

where  A  =  A(k, u)  is  to  be  determined,  and  7  =  (c 02/co2  —  k2)^.  The  radiation 
condition  is  satisfied  by  requiring  the  square  root  (cu2/c02  —  k2)^  to  have  the  same 
sign  as  u  when  the  argument  is  positive  and  to  be  positive  imaginary  when  the 
argument  is  negative. 

Within  the  bubble  layer  it  is  assumed  that  pressure  fluctuations  satisfy  the 
wave  equation  (6.2.1)  with  cQ  replaced  by  c,  so  that  we  can  write 

p  =  (fie^3  +  C'e-^3)ei(k  x-^,  7  =  (w2/c2  -  k2)  3,  (6.4.5) 

where  the  B(k,uj),  C(k,  u)  are  to  be  determined. 

The  coefficients  A,  B,  C  are  obtained  by  application  of  conditions  of  continuity 
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of  pressure  at  x3  =  0,  h  and  (when  the  mean  flow  Mach  number  is  infinitesimal) 
continuity  of  the  x3  component  of  velocity.  The  pressure  conditions  yield 


F(u)  =  B  +  C,  at  x3  =  0, 

A  =  Be^h  +  Ce~m,  at  X3  =  h. 


(6.4.6) 


Continuity  of  normal  velocity  requires  that  (1/ p)dp/dx3  be  continuous  at 
x3  =  h,  i.e.  that 

^-A  =  Bem-Ce~m.  (6.4.7) 

7  Po 

The  solution  of  equations  (6.4.6),  (6.4.7)  give 


A 


F(u) 


cos(y h)  -i  (£)sin(7 h) 


(6.4.8) 


The  substitution  of  this  result  into  equation  (6.4.4)  and  application  of  the  integral 
operator  (6.4.3)  the  supplies  the  acoustic  pressure  in  the  water  in  the  form  [43] 


P(x,t) 


1  r  r°°  f(u)  e4k-x+7C3 -h.)-ut)  d2kdu 

(2 7r)2  J J-00  cos(7 h)  —  i  sin(7h) 


(6.4.9) 


6.4.3  The  far  field  acoustic  pressure 

The  radiated  acoustic  power  is  determined  by  the  far  field  limit  of  the  pressure 
determined  by  equation  (6.4.13).  In  the  propagating  region  at  x  (|x|  — *  00) 
the  integral  is  dominated  by  contributions  from  those  wavenumbers  k  in  the 
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neighborhood  of  (ki,  £+7)  —  /c0x/|x|,  and  the  leading  approximation  to  the  integral 
(the  stationary  phase  approximation  [12,40,42,44])  is  obtained  by  expanding  the 
integrand  about  this  point.  This  gives 


P(x,t)  ^ 


cos  9 


27TcJx 


—iujF(uj)  e  |x| /Co)du  x3  ,  ,  (a  a 

-,  cosd  =  1 — -,  |x|  — >•  00,  (6.4.10) 


o|xl  J  00  cos(7  h)  —  i  (  )  sin  (7/1) 


where  the  denominator  of  the  integrand  is  evaluated  at  (k\ ,  ^2,7)  =  /c0x/|x|.  Thus, 
7  =  ka  cos  9  and 

7  =  ~f(u,9), 

Co 


where 


/M)  = 


0  I  2  nii  ■ 

—  COS  0+11 - 2  —  lrl  — 


wt 


-1 


1/2 


(6.4.11) 


Jb 

The  term  (c0/c0)2  cos2  6  in  this  formula  can  be  neglected  except  at  large  frequencies 
when  uj  is  comparable  to  c +,(c0/c0).  At  very  low  frequencies  /  ~  1. 


Hence,  equation  (6.4.10)  reduces  to 


P(x,t)  ^ 


cos  6  r°°  —  iuF(cu)e  du 

27TC0|x|  J Loo  cos  (^)  -  i{Z / f)  cos  9 sin  (^)  ’ 


X  — >•  oo. 


(6.4.12) 


where  /  =  f(u>,9). 

As  h  — >  0  the  integral  representation  (6.4.12)  reduces  to  the  exact  result  (6.2.2) 
in  the  absence  of  the  bubble  layer;  similarly  both  of  the  approximate  representations 
(6.4.10),  (6.4.12)  reduce  to  the  corresponding  far  field  approximation  (6.2.3).  Note 
also  that  whereas  the  radiation  in  the  absence  of  the  bubble  layer  has  the  usual 
dipole  directivity,  the  peak  sound  radiating  in  the  direction  normal  to  the  interface, 
the  corresponding  directivity  in  the  presence  of  the  bubble  layer  is  frequency 
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dependent.  In  particular  at  those  frequencies  at  which  cos (uh/c)  =  0  it  is  clear 
that  the  radiation  directivity  becomes  omnidirectional.  This  is  because  at  these 
frequencies  the  collective  response  of  the  bubbles  in  the  vicinity  of  the  jet  impact 
point  produces  near  resonant  volume  fluctuations  of  the  layer,  which  accordingly 
radiates  as  a  monopole  source. 


6.4.4  The  acoustic  pressure  frequency  spectrum 


The  definitions  (6.2.4),  (6.2.5)  permit  the  acoustic  pressure  frequency  spectrum  to 
be  cast  in  the  form 


4>pP(x,a;) 


cos2  6  k°&FF(u) 


4vr2|x|2  [cos(f  /)  -i(Z/)cos0sin(f  /)] 


2  ’ 


x  — y  oo. 


(6.4.13) 


This  reduces  to  Foley’s  [1,2]  formula  (6.2.6)  when  the  nominal  thickness  of  the 
bubbly  interface  h  — y  0.  The  non-dimensional  Foley  spectrum 


4>Pp(x,q;)t/Z>3  _  cos 2  9  a(k0D )2  16  4  141 

FKD/ |x|)2  ~  2vr2  (1  +  e*(uD/UY)l  K  '  '  ] 

is  plotted  as  the  broken  line  curve  ( - )  in  Figure  6.22  for  radiation  in  the  normal 

direction  to  the  interface  (6  =  0)  in  the  case 


U  =  27.1m/s,  D  —  18.6mm,  a—  17.28.  (6.4.15) 


The  corresponding  spectra  in  the  presence  of  a  bubble  layer  of  thickness  h  =  10 
mm  are  plotted  as  the  solid  line  curves  on  Figure  6.22  for  cases  where  the  mean 
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bubble  resonance  frequency  ft,  =  ojb/ 2tt  =  5, 10,  20,  30, 40  kHz,  and  for 

ca  =  100 m/s,  p  =  980 kg/m3.  (6.4.16) 


All  of  these  spectra  exhibit  a  large  peak  at  k0D  0.2,  or  at  about  2.5  kHz  (typically 
10-20  dB  above  the  Foley  spectrum).  This  is  the  fundamental  resonance  frequency 
of  the  ’quarter  wave’  series 


i oh  unh 
Co  C0 


n  =  0,1,2,..., 


(6.4.17) 


at  which  cos (coh/c0)  =  0  for  a  tube  of  open  length  h  open  at  one  end  (the  upper 
interface  at  x 3  =  h  corresponds  to  the  ’closed’  end).  The  only  relevant  additional 
peaks  in  Figure  6.22  correspond  to  the  resonances  n  —  1  and  2.  It  is  also  clear  from 
the  figure  that  propagation  through  the  bubble  layer  is  progressively  attenuated  at 
frequencies  approaching  and  exceeding  the  mean  bubble  resonance  frequency  fb- 
This  eventually  cuts  off  the  higher  frequency  components  of  the  sound  that  radiate 
in  the  absence  of  bubbles. 


Figure  6.22:  The  non-dimensional  acoustic  pressure  frequency  spectrum.  The  dashed 
line  is  for  h  =  0  and  the  solid  line  is  for  a  bubble  layer  of  thickness  0.39". 
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This  very  large  amplification  of  the  radiated  sound  by  the  bubble  layer  is 
probably  unrealistic,  because  it  is  unlikely  that  the  layer  will  be  stable  enough  to 
maintain  the  constant  value  of  h  required  for  a  resonant  mode  to  saturate.  It  is 
only  the  relatively  slowly  varying  components  of  h  that  are  likely  to  contribute 
significantly  to  a  resonance-like  behavior.  The  magnitude  of  such  terms  can  be 
estimated  by  averaging  the  spectrum  (6.4.13)  with  respect  to  its  dependence  on  h. 

The  simplest  hypothesis  is  to  assume  that  h  is  uniformly  distributed  in 
probability  between,  say,  zero  and  a  maximum  value  fimax.  The  rectified  form  of  the 
acoustic  pressure  frequency  spectrum  (6.4.13)  is  then  obtained  by  integrating  with 
respect  to  h  over  (0,  hmax)  and  dividing  by  hmax.  The  result  can  be  written  as 


$pP(x,o;) 


cos  9  2  A  (Q,9) 


47T2  |x  | 


A 


|x|  — »  oo,  Id  = 


Ujllr 


(6.4.18) 


where 


A(n,6) 


Z  cos  9  fhuiax 


dh 


i(Z/f)  cos  9  sin 


(6.4.19) 


The  factor  Z  cos  9  before  the  integral  is  included  because  the  value  of  the  integral 
tends  to  be  dominated  by  contributions  from  the  peaks  of  the  integrand  (seen  in 
Figure  6.23)  causing  its  order  of  magnitude  1  /{Z/  cosd).  Figure  6.23  shows  the 
plot  of  A(f2,  9)  as  a  function  of  11  =  cohmax/c0  for  fimax  =  10  mm,  9  =  0°,  and  for  the 
mean  bubble  resonance  frequencies  fb  =  5, 10,  20,  30,40  kHz. 

The  function  A(£2,9)  turns  out  to  be  only  weakly  dependent  on  the  radiation 
direction  9  when  1  <  <  10,  and  peaks  roughly  at  ojhmax/c0  ~  7t/2.  In  this  range 

the  directivity  of  the  rectified  acoustic  pressure  frequency  spectrum  (6.4.18)  is 


109 


approximately  proportional  to  cos  9,  which  is  midway  between  the  dipole  obtained 
by  Foley  et  al  [7,8]  in  the  absence  of  a  bubble  layer  and  the  omnidirectional  monopole 
that  is  predicted  by  Equation  (6.4.13)  at  this  quarter-wave  resonance  frequency 
for  a  constant,  uniform  bubble  layer.  At  lower  frequencies  A (£2,9)  ~  Zcos9  and 
Equation  (6.4.18)  provides  a  smooth  transition  from  the  cos  9  to  cos 2  9  dependence 
of  the  directivity. 


Figure  6.23:  Illustrating  the  variation  of  A(Q,  9)  when  9  =  0°,  p  =  980kg/m3, 
c0  =  100  m/s,  hmax  =  10  mm  and  fb  =  5, 10,  20,  30, 40  kHz. 


Figure  6.24  compares  the  Foley  [7,8]  spectrum  and  the  rectified  spectrum  for 
radiation  in  the  normal  direction  (9  =  0)  for  conditions  (6. 4. 13), (6. 4. 18)  when 
/imax  =  10  mm  and  for  bubble  resonance  frequencies  of  fb  =  5, 10,20,30,40  kHz. 

In  all  cases  there  are  significant  increases  in  the  sound  pressure  level  in  the  range 
0.1  <  k0D  <  1  as  large  as  10  dB  when  fb  exceeds  20  kHz.  At  higher  frequencies  the 
bubble  layer  causes  the  levels  to  fall  rapidly  below  that  in  the  absence  of  bubbles. 
For  the  range  of  frequencies  illustrated  in  Figure  6.24  the  high  frequency  decay  is 
<FpP(x,  uf)  «  (k0D)~3/2. 
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Figure  6.24:  The  non-dimensional  acoustic  pressure  frequency  spectrum  in  the  nor¬ 
mal  direction  6  =  0  and  conditions  (6. 4. 16), (6. 4. 18).  The  dotted  line  represents 
Foley’s  spectrum  [7,8]  for  h  =  0.  The  solid  lines  are  of  the  rectified  acoustic  spec¬ 
trum  for  hmax  =  10  mm  and  ft,  =  5, 10, 20,  30, 40  kHz. 


6.5  Conclusions 

The  motivation  behind  these  mathematical  models  came  from  observations  of 
bubble  entrainment  due  to  gas  jets  impacting  a  free  air-water  interface.  The 
resulting  model  predicts  that  the  radiated  sound  will  exhibit  a  dipole  radiation 
pattern,  with  or  without  a  bubbles  at  the  interface.  At  a  range  of  frequencies, 
the  acoustic  spectrum  is  increased  by  10-30  dB.  At  higher  frequencies  which  have 
wavelengths  that  are  on  the  same  length  scale  of  the  mean  cloud  radius  R  or  bubble 
layer  thickness  h,  the  sound  power  is  predicted  to  be  lower  than  in  the  absence  of 
bubbles. 

The  sum  of  the  measured  unsteady  force  spectrums  from  all  the  force  sensors 
were  used  as  an  input  to  Foley’s  acoustic  model.  The  predicted  levels  were  generally 
about  15  dB  below  the  measured  acoustic  spectrum.  This  indicates  that  the 
current  acoustic  model  is  not  representing  all  of  the  energy  transmitted  through 
the  air-water  interface.  Some  of  this  could  be  in-air  jet  noise  transmitting  directly 
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through  the  interface.  The  inaccuracy  in  the  model  could  also  stem  from  the 
assumption  that  the  jet  impact  area  can  be  idealized  as  a  delta  function.  In  reality, 
the  jet  impacts  the  surface  over  a  finite  area  that  could  lead  to  contributing  sources 
over  the  whole  impact  area. 

Finally,  the  sum  of  the  measured  unsteady  force  spectrnms  from  all  the  force 
sensors  were  used  as  an  input  to  the  bubble  hemisphere  model  The  predicted  levels 
for  the  1/16"  nozzle  was  about  20  dB  below  the  measured  levels  for  most  of  the 
spectrum.  While  the  overall  levels  are  lower,  the  shape  of  the  model  matches  the 
measured  spectrum.  The  predicted  levels  for  the  1/8"  nozzle  was  about  10  dB 
below  the  measured  levels  over  most  of  the  spectrum.  Again,  the  overall  predicted 
levels  were  lower,  but  the  shape  of  the  model  matches  the  measured  data.  There 
are  a  few  possible  explanations  for  this  discrepancy.  The  in-air  jet  noise  transmitted 
directly  through  the  interface  is  not  accounted  for  in  the  model.  Also  certain 
parameters,  such  as  the  void  fraction  and  bubble  size  within  the  bubble  hemisphere 
were  estimated.  This  would  effect  the  sound  speed  within  the  bubble  hemisphere, 
and  would  change  the  predicted  level  of  the  model.  Direct  measurement  of  the 
composition  of  the  bubble  cloud,  including  the  void  fraction  and  bubble  size,  could 
improve  the  quality  of  the  bubble  model. 


Chapter  7 


Conclusions 

This  chapter  acts  as  a  link  between  the  work  performed  in  this  investigation  and 
future  research.  Section  7.1  summarizes  the  contents  of  this  thesis  along  with  the 
most  important  results.  Section  7.2  provides  recommendations  to  vehicle  designers 
for  reducing  the  self  noise  of  supercavitating  vehicles.  Finally,  Section  7.3  discusses 
the  questions  brought  up  by  this  study,  and  proposes  some  future  experiments. 

7.1  Summary 

This  thesis  describes  the  investigation  into  the  self  noise  from  a  ventilated  super¬ 
cavitating  vehicle.  This  problem  was  approached  analytically  and  experimentally. 
Chapter  2  extended  Foley’s  original  theory  and  derived  acoustic  models  for  sound 
traveling  through  the  air-water  interface  in  the  presence  of  either  a  bubbly  layer  or 
hemisphere.  Chapter  3  discussed  measurements  of  the  gas  jet  system  to  be  used  in 
later  experiments,  ft  was  important  to  understand  the  flow  properties  and  radiated 
noise  of  the  jet  system.  This  turned  out  to  be  the  major  sources  of  uncertainty 
in  later  experiments.  Chapter  4  described  the  measurements  of  the  unsteady 
force  spectrum  of  a  gas  jet  hitting  a  rigid  wall.  These  results  were  compared  to 
measurements  taken  by  Foley.  The  spectrums  measured  in  this  experiment  were 
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similar  to  those  measured  by  Foley,  but  with  less  noise.  Chapter  5  outlined  the 
acoustic  measurements  of  gas  jets  hitting  the  air-water  interface  with  no  bubble 
entrainment.  This  experiment  matches  with  the  original  acoustic  model  developed 
by  Foley.  The  measured  force  data  was  used  as  an  input  to  Foley’s  model,  but  the 
predicted  levels  were  much  lower  than  those  measured.  The  relationship  between 
jet  exit  velocity  and  acoustic  pressure  levels  was  also  explored  in  this  chapter.  This 
led  to  the  observation  that  the  acoustic  pressure  level  is  also  dependent  on  the 
jet  impact  area  at  the  surface.  Gas  jets  with  larger  impact  areas  on  the  air-water 
interface  had  a  higher  slope  power  law.  Chapter  6  described  the  acoustic  pressure 
spectrum  measurements  for  the  bubble  entrainment  cases.  Observed  dimensions  of 
the  entrained  bubble  clouds  were  used  as  inputs  to  the  acoustic  model  developed  in 
Chapter  2.  The  acoustic  model  predicted  an  amplification  in  acoustic  signal  over  a 
range  of  frequencies  that  was  verified  by  the  measured  data. 


7.2  Recommendations  to  reduce  self  noise 

The  most  important  factor  in  the  self  noise  due  to  the  ventilation  system  of 
supercavitiating  vehicles  is  the  exit  velocity  of  the  gas  jets.  Higher  velocity  jets  lead 
to  higher  radiated  noise  (Figures  4.4-4. 6,  Section  4.3.1).  A  relationship  between  the 
size  of  the  jet  impact  area  and  radiated  noise  was  also  discovered.  Larger  diameter 
jet  impact  areas  are  more  sensitive  to  changes  in  exit  velocity  (Section  4.3.1).  It 
seems  that  the  best  solution  would  be  to  vent  the  minimum  amount  of  gas  to 
produce  a  stable  cavity  around  the  vehicles.  It  was  shown  that  in-air  tones  from 
the  ventilation  system  could  be  transmitted  through  the  air-water  interface  (Section 
4.2.2).  Detailed  in-air  measurements  of  the  ventilation  system  would  be  advisable. 
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It  would  be  a  good  idea  to  keep  the  ventilation  flow  as  smooth  as  possible  to 
eliminate  this  self  noise  (Section  2.3.2).  It  has  also  been  shown  that  the  existence  of 
bubbles  can  increase  an  acoustic  signal  (Equations  6.3.15  and  6.3.16,  Figures  6.16, 
6.19,  6.20).  Therefore  any  ventilation  strategies  that  visibly  produces  bubbles  at  the 
interface  should  be  avoided  if  possible.  But  in  a  realistic  saltwater  environment,  the 
complete  elimination  of  bubbles  near  a  high  speed  underwater  vehicle  is  impossible. 

7.3  Recommendations  of  future  work 

The  main  recommendation  from  this  thesis  would  be  to  study  the  interactions  of 
gas  jets  with  deflector  plates  above  the  air-water  interface.  This  topic  was  just 
beyond  the  scope  of  this  investigation  to  fully  cover.  But  interesting  analytical 
predictions  about  the  influence  of  a  deflector  were  derived.  Begin  by  writing  the 
outgoing  acoustic  pressure  in  terms  of  the  Fourier  integral 

f‘  COO 

p(x,f)=/  /  p(k,cu)ei(k'x+7*3-“Vkdu  (7.3.1) 

J  J — oo 

where  k  =  (ki,k2,  0)  and  p(k,o;)  is  the  space-time  Fourier  transform  of  the  pressure 
perturbation  on  the  interface  X3  =  0.  This  can  be  rearranged  to  form 

p(k,co)  =  F(t)5(x1)5(x2)e-^-^dx1dx2dt  =  ^  (7.3.2) 

(27 r)J  J  (2vr)^ 

with  7  =  {oj2  /  c2  —  k2)^ .  The  radiation  condition  is  satisfied  by  requiring  (c o2  /  c2  —  k2)^ 
to  have  the  same  sign  as  to  when  u2 /c20  >  k 2  and  to  be  positive  imaginary  for  other 
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real  values  of  k.  Therefore 


p(x,t)  =  — [  [  F(u)ei^+^-“t)d2\Ldu,Xz>U.  (7.3.3) 

In  the  far  field  at  x(|x|  — y  oo )  the  wavenumber  integral  is  dominated  by  the  region 
near  (ki,  &2, 7)  =  /c0x/|x|,  and  the  leading  approximation  to  the  integral  (the 
’stationary  phase  approximation’)  is  obtained  by  expanding  the  integrand  about 
this  point,  yielding 


eosd  r°° 

p(x,t)  «  — -Vr  /  -icuF(cu)e-^(*-|x|/Co)dcu,  |x|  -)■  00  (7.3.4) 

27TC0|x|  J _ oo 

This  implies  that  p(x,cu)  «  — icaF(a;)e*u;lx^Co  cos 9/{2'kcg |x|)  when  |x|  — >  00,  which 
leads  back  to  Foley’s  original  theory  [7,8]. 


Figure  7.1:  Diagram  of  gas  jet  interacting  with  a  deflector  plate  normal  to  the 
interface. 

Next  consider  the  case  in  which  a  rigid  deflector  plate  has  been  placed  between 
the  jet  and  the  water  at  distance  h  from  the  surface,  as  shown  in  Figure  7.1.  In 


a  first  approximation  it  may  be  assumed  that  the  unsteady  force  experienced  at 
the  plate  is  the  same  as  that  F(t)  produced  by  impingement  on  the  gas-water 
interface.  The  impact  with  the  plate  generates  dipole  aerodynamic  sound.  The 
direct  radiation  pd  within  the  cavity  is  given  by  Curie’s  formula 


Pd 


d  r  1 

dx'3  47r|x' 


F 


(7.3.5) 


where  x'  =  (. Xi,x2,x3  —  h )  and  C\  is  the  speed  of  sound  within  the  cavity.  In  the 
region  0  <  x':i  <  h  between  the  deflector  and  the  interface  the  result  (7.2.5)  can  also 
be  expressed  in  the  following  manner  corresponding  to  equation  (7.2.3): 


1 

(2vr)2 


-F(co)ei(k  x'+yx^-ut)d2kdu,x'3  >  0, 


(7.3.6) 


where  y'  =  (cu2 / c2  —  k2)^ .  The  incidence  on  the  interface  £3  =  h  of  the  acoustic 
pressure  pd  produces  (to  first  order)  a  reflected  wave  of  equal  amplitude  and  phase. 
The  Fourier  transform  of  the  net  interface  pressure  is  therefore 


P{  k,w) 


F{pj)ell'h 

(2vr)2 


(7.3.7) 


The  analogues  of  equations  (7.2.2)  and  (7.2.3)  accordingly  yield  the  acoustic 
pressure  in  the  water  in  the  form 


F(u)e 


P(x,t) 


1 


OO 


,i  (k- x+7^3  +7'  h— cot )  d2  kdcd 


(7.3.8) 
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which  can  be  approximated  as 


P(x,t) 


so  that 


cos  0  r°° 

2vrc0|x|  J ^ 


-zo;r  (cu)e  vci  c°  /  du, 


p(x,w) 


-iujF(uj)  COS  6  i<jj\-x\/c0-\-iwh( \ - 2 

v  -e  v  ci  c° 


2  2 


27tcJx| 


x|  — *  oo,  X3  >  0 

(7.3.9) 

(7.3.10) 


Because  4  >  it  follows  that  the  acoustic  pressure  is  given  by 


^>»®(x,  u>) 


cos2  6 
47T2  |x  | 2 


kl&FF^tu),  |x|  — *  OO 


(7.3.11) 


This  is  the  same  result  as  if  the  jet  were  to  be  impinging  directly  on  the  interface. 

Next  consider  a  deflector  plate  oriented  within  a  cavity  that  is  oriented  so  that 
its  outward  normal  is  in  the  forward  direction,  as  shown  in  Figure  7.2. 


Figure  7.2:  Diagram  of  gas  jet  interacting  with  a  deflector  in  a  cavity  oriented  in  the 
forward  direction. 
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The  dipole  source  at  the  deflector  is  parallel  to  the  X2  direction  of  Figure  7.1, 
and  equations  (7.2.9)  and  (7.2.10)  become 


P(x,t) 


1 

w? 


Y(k) 


(7.3.12) 


p(x,t) 


x2x3 


27rc0|x|3y/ (c0/ci)2  —  sin2  8  J~c 


-lur  (u)e  ei 


du,  |x|  — *  00,  x3  >  0. 
(7.3.13) 


Therefore 

$pp(x,o;)  «  ,  —  kl$FF(u),  | x |  ->■  00  (7.3.14) 

47r2|x|6(|f  —  sm2  9) 

The  sound  speed  ratio  c0/c\  ~  5.  The  acoustic  pressure  spectrum  in  the  water 
is  therefore  reduced  by  about  10  x  log10(25)  ~  14  dB  compared  to  a  deflector 
oriented  out  to  the  interface.  Taking  in  air  measurements  of  gas  jets  impinging  on 
deflector  plates  at  various  angles  would  be  a  good  start  towards  understanding  this 
interaction.  The  next  step  would  be  to  take  underwater  measurements  of  gas  jets 
hitting  a  deflector  plate  above  the  interface. 

Improved  measurements  of  the  composition  of  the  bubble  hemisphere  would 
likely  improve  the  quality  of  the  bubble  hemisphere  model.  This  would  include  direct 
measurements  of  the  void  fraction  and  bubble  size  within  the  bubble  hemisphere.  It 
would  also  be  interesting  to  study  the  indentations  of  the  low  speed  gas  jets  hitting 
the  air-water  interface.  In  these  experiments,  the  water  was  still.  If  the  water  were 
moving  at  a  steady  rate,  the  indentation  might  shear  off  and  entrain  bubbles.  This 
could  potentially  produce  a  layer  of  bubbles  near  the  interface.  It  could  be  useful 
to  determine  the  minimum  indentation  size  necessary  to  entrain  bubbles  in  moving 


water. 


Appendix  A 


Review  of  select  topics  in  linear  acoustics 

This  appendix  summarizes  some  of  the  fundamentals  of  acoustic  theory  related  to 
this  dissertation.  Material  is  taken  from  Blackstock[45],  Howe  [12],  and  Pierce [44], 
See  these  references  for  greater  detail  on  each  topic  or  a  more  general  review  of  the 
principles  of  linear  acoustics. 


A.l  The  linear  wave  equation 


The  acoustic  wage  equation  is  derived  from  the  conservation  and  constitutive 
equations  which  describe  a  fluid.  We  consider  a  fluid  volume  V  defined  by  a  surface 
S  in  which  body  forces  and  thermal  and  viscous  effects  are  negligible.  The  fluid 
inside  the  volume  has  velocity  u(x,  t),  pressure  P(x,  t),  and  density  p(x,  t)  at  a 
position  x  =  (aq,  aq,  X3)  and  time  t. 

Conservation  of  mass  requires  that  the  increase  of  mass  inside  V  must  be  equal 
to  the  influx  of  mass  through  S.  This  is  expressed  by  the  continuity  equation 


Dp 

Dt 


+  pV  •  u 
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where  the  material  derivative  D/Dt  —  d/dt  +  u  •  V. 

Conservation  of  momentum  holds  that  the  time  rate  of  change  of  momentum 
inside  V  is  equal  to  the  sum  of  all  forces  acting  on  the  fluid  volume  and  the 
momentum  flux  through  the  surface  S.  Neglecting  body  and  viscous  forces, 
momentum  conservation  can  be  expressed  by 


Du 

pm+VP 


o 


Finally  we  require  an  equation  of  state,  a  relationship  between  thermodynamic 
variables.  We  will  assume  that  our  fluid  volume  is  isentropic,  the  entropy  s  in  the 
volume  is  constant. 

p  =  p(P) 

We  consider  that  the  pressure  P  and  density  p  can  be  decomposed  into  a  mean,  or 
static,  and  excess  constituent  such  that 


P  =  Po+P 


p  =  p0  + dp 


where  pQ  and  pQ  are  respectively  the  mean  pressure  and  density  in  the  fluid  volume 
5p  is  the  excess  density,  and  p  is  the  acoustic  pressure.  This  leads  to  the  Taylor 
expansion  of  the  equation  of  state 


Pc 


+ 


p 


Pi 
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where  A,  B ,  etc.  are  empirically  determined  coefficients.  Introduce  the  sound  speed 
c  in  the  fluid  defined  by 


<,  dP  A  B  p  —  p0 

c2  =  —  =  —  +  —  - - —  +  •  •  • 

dp  po  po  po 

Letting  p  — *  p0,  we  find  that  c  becomes  a  constant  cQ  and  that  A  =  pDc2.  Substitute 
for  p  and  p  to  get  the  equation  of  state 


P  =  c205p  1  + 


B  5p 
VAjo 


+  ■■■ 


To  simplify  these  equations,  we  first  observe  that  the  changes  in  the  fluid 
properties  due  to  sound  energy  are  small.  In  a  quiet  fluid,  P  =  pQ,  p  =  pa,  and  u  =  0 
across  V.  We  also  assume  that  the  fluid  is  homogeneous  with  constant  pressure  and 
density  across  the  volume.  When  sound  is  introduced  the  excess  pressure,  density, 
and  induced  particle  velocity  are  very  small  compared  to  their  static  values 


Sp  <C  p 

\P\  <  PoC2 
\u\  <C  Co 

Substitute  for  P  and  p  and  expand  the  continuity  equation  as 


d_ 

dt 


(Sp)  +  uWSp  +  poV  ■  u  +  SpW  •  u  =  0 


Considering  our  small  value  approximations,  we  note  that  the  second  and  fourth 
terms  on  the  left  hand  side  above  are  an  order  of  magnitude  smaller  than  the  first 
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and  third.  Eliminating  them  leads  to  the  linearized  continuity  equation 


+PoV  •  u  =  0 


Next  apply  the  same  procedure  to  the  conservation  of  momentum  equation. 


<9u  <9u 

Povw  +  PoU  •  Vu  +  dp—  +  dp u  ■  Vu  +  Vp  =  0 

at  at 


As  before,  keep  only  the  first  and  third  terms  of  the  series.  The  linearized 
conservation  of  momentum  equation  is 


<9u 

P°m+Vp  =  0 


The  equation  of  state  is  linearized  to 


P  =  c20dp 


Use  this  to  eliminate  dp  from  the  continuity  equation 


Id p 

"T  +  PoV  •  u  =  0 
ci  at 


Taking  the  difference  of  the  time  derivative  of  this  equation  with  the  spatial 
derivative  of  the  linearized  momentum  equation  results  in  the  linear  acoustic  wave 
equation 


\Cq  d t2 


p  =  0 


It  is  often  convenient  to  consider  sound  transmission  in  the  frequency  domain  rather 
than  the  time  domain.  The  wave  equation  is  brought  into  the  frequency  spectrum 
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by  applying  a  Fourier  transform,  which  gives  the  Helmholtz  equation 

(kl  +  v2)p  =  0 

where  the  acoustic  wavenumber  ka  =  oj/c0  and  u  is  angular  frequency.  Next 
consider  a  fluid  volume  that  contains  an  acoustic  source.  The  source  can  be 
considered  to  add  a  volume  source  to  the  right  side  of  the  continuity  equation  and 
a  force  term  to  the  right  side  of  the  momentum  equation.  Carrying  these  terms 
through  the  derivation  yields 

(^-V2)p  =  fXx,i) 

where  the  volume  and  force  terms  have  been  replaced  by  a  generalized  pressure 
source  F(x,t).  This  result  is  called  the  inhomogeneous  wave  equation. 

In  general,  the  solution  of  the  wave  equation  and  Helmholtz  equation  can  be 
decomposed  into  an  outgoing  and  incoming  wave.  At  time  t  —  0  an  acoustic  source 
begins  to  generate  sound  in  the  fluid.  Causality  requires  that  for  t  >  0  the  wave 
energy  radiates  outward,  and  in  the  absence  of  other  sources  or  bodies,  there  is  no 
incoming  energy.  The  consideration  that  all  incoming  solutions  of  the  wave  equation 
must  be  zero  is  called  the  radiation  condition. 


A. 2  Acoustic  impedance,  intensity,  and  power 

Acoustic  impedance  quantifies  the  sound  absorption  of  a  fluid,  It  is  defined  as  the 
ratio  of  acoustic  pressure  to  particle  velocity  in  the  medium.  For  an  outgoing  wave, 
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the  characteristic  impedance  ZQ  is  given  as 


v  _P  _ 

Z o  PoCo 

u 


A  medium  is  considered  acoustically  rigid  if  its  impedance  Z  — >■  oo.  A  normally 
incident  plane  wave  traveling  through  a  medium  of  impedance  Z\  is  reflected  back 
entirely  when  it  reaches  a  medium  of  impedance  A2  when  Z\  <C  Z2.  If  Z\  S>  Z2, 
there  is  again  no  transmission.  The  reflective  wave  identically  cancels  the  incident 
wave  and  the  interface  is  known  as  a  pressure  release  surface,  because  the  acoustic 
pressure  p  is  zero  at  the  interface.  When  the  impedance  Z  — >  0  the  medium  is 
considered  acoustically  compliant. 

Acoustic  intensity  is  the  time  average  of  sound  energy  flow  through  a  unit  area 
in  the  direction  normal  to  that  area.  In  integral  form,  acoustic  intensity  /  is 


r^av 

I  =  —  /  pudt 

tav  JO 


The  directivity  of  radiated  sound  is  the  directional  dependence  of  the  acoustic 
intensity.  For  sound  radiating  omnidirectionally,  intensity  on  a  spherical  surface  of 
radius  r  simplifies  to  a  scalar  quantity 


/  = 


p(r): 

PoCo 


The  acoustic  power  W  passing  through  a  surface  S  is  the  integral  of  the  sound 


W  =  /  IdS' 

Js 


intensity  over  that  surface 
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where  S'  is  a  differential  element  of  the  surface  S.  For  omnidirectional  sound,  the 
power  through  a  sphere  of  radius  r  is 

W  =  4t rr2/ 

A. 3  Acoustic  compactness  and  the  near  and  far  fields 

Am  acoustic  source  occupying  a  fluid  region  of  characteristic  dimension  l  and 
angular  frequency  oj  radiates  sound  of  wavelength  2tt cq/uj.  The  source  is  considered 
to  be  acoustically  compact  when  l  is  much  smaller  than  this  wavelength,  or 
u)l / cQ  =  kQl  <£.  1.  The  acoustic  far  held  exists  at  point  |x|  many  wavelengths  from 
the  source  where  k0 |x|  3>  1.  The  acoustic  near  held  is  the  region  in  the  immediate 
vicinity  of  the  source  where  kQ  <  1. 

A. 4  Monopoles,  dipoles,  and  quadrupoles 

A  monopole  is  an  acoustic  source  which  can  be  considered  analogous  to  a  pulsating 
sphere.  Sound  radiates  from  the  source  omnidirectionally,  diminishing  in  the  far 
held  (as  |x|  — y  oo)  like  l/|x|.  A  general  expression  for  a  monopole  is  a  source  of  the 
form  q(t)S(x),  where  q(t)  is  the  strength  of  the  source. 

A  dipole  is  a  source  of  the  form 

AywsM] 

with  strength  fj(t).  A  dipole  is  analogous  to  an  oscillating  sphere.  Sound 
emanates  from  the  source  as  two  lobes  protruding  in  the  direction  of  that  os- 
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dilation  as  shown  in  Figure  A.l,  diminishing  as  1  /|x|2  in  the  far  held.  A  point 
dipole  is  also  equivalent  to  two  monopoles  of  opposite  strength  a  short  distance  apart. 


Dipole  Directivity  Quadrupole 

Directivity 


Figure  A.l:  Directivity  of  the  sound  from  a  dipole  and  quadrupole  source. 

A  source  distribution  equivalent  to  four  monopoles  with  net  strength  zero  is 
called  a  quadrupole.  A  quadrupole  source  is  expressed  as 

d2T 

^■t:> 

and  diminishes  in  the  far  held  as  1  / 1 x | 3 .  The  directivity  of  a  quadrupole  is  shown 
in  Figure  A.l. 

An  acoustic  source  can  be  represented  as  the  sum  of  component  sources  of 
monopole,  dipole,  and  quadrupole  strength.  Typically  it  is  necessary  only  to 
consider  the  leading  order  term  in  determining  the  radiated  acoustic  pressure 
generated  by  the  source. 


Appendix  B 


Acoustic  Green’s  function 


This  appendix  discusses  the  use  of  acoustic  Green’s  functions  as  solutions  to  the 
inhomogeneous  wave  equation.  The  material  is  taken  from  Howe’s  Acoustics  of 
Fluid  Structure  Interactions^]  and  Theory  of  Vortex  Sound] 46]. 

The  free  space  Green’s  function  is  the  outgoing  solution  of  the  wave  equation 
generated  by  the  point  source  <5(x  —  y )5{t  —  r)  in  an  unbounded  fluid  located  at  the 
point  x  =  y  and  at  time  t  —  r. 

(^-V2)G  =  S(x-y)6(t-T) 
where  G  =  0  for  t  <  r.  This  has  the  solution 


G(x,y,t-r) 


1 

47t|x  —  y 


r 


x-yh 

<2  ) 


The  inhomogeneous  wave  equation 


\  c20  d t2 


V  =  F(x,t) 
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describes  the  behavior  of  waves  generated  by  a  distributed  source  F(x,  t).  This 
source  can  be  regarded  as  an  array  of  point  sources,  such  that 


F(x,f)  = 


F{ y,  t)S(x  -  y)S(t  -  r)d3ydr 


The  solution  for  each  constituent  source  of  strength  F( y,  t)S(x  —  y )5(t  —  r)d3ydr 
is  F(y,r)G(x,y,t  —  r)<f3ydr.  Add  these  individual  contributions  to  obtain 


p(x,f)  = 


p(x,t)  = 


47 T 


F{ y,  t)G(x,  y,  t  -  r)d3ydr 


F(y,r,slt_T_i2^y1]d3ydT 


F(y,t 


lx — y I ' 


-d3y 


^  J- oo  lx  —  y  I 

Green’s  functions  can  be  applied  to  more  general  acoustic  problems  involving 
solid  or  fluid  boundaries,  and  have  equivalent  expressions  in  the  time  and  frequency 
domains.  The  use  of  the  Green’s  function  will  be  demonstrated  by  solving  the 
inhomogeneous  Helmholtz  equation.  Consider  an  acoustic  source  on  the  boundary 
of  a  volume  of  fluid  defined  by  some  surface  S.  The  source  generates  acoustic 
pressure  p(x,  oo)  within  the  fluid  volume  which  satisfies 


(V2  +  kl)p  =  0 


subject  to  the  boundary  condition  p  =  p'(x,u)  on  S.  Introduce  a  Green’s  function 
G(x,lo)  with  outgoing  wave  behavior 


(V2  +  k20)G  =  0 
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with  the  boundary  condition  G 


G'(x,u)  on  S.  Taking  the  difference  of  G  and  p 


pV2G  -  GV2p  =  0 


Integrating  over  the  fluid  volume  V, 


(pV2G  -  GV2p)dV 


0 


Introduce  the  vector  identities 


pV2G  =  V2(PG)  -  V  •  (GVp)  -  VG  •  Vp 
GV2p  =  W2{Gp)  -  V  •  (pVG)  -  Vp  •  VG 


into  the  integrand 

[  [V  •  (pVG)  -  V  ■  (GVP)}dV  =  0 

Jv 

Applying  the  divergence  theorem  reduces  this  to 


[V  •  (pVG)  -  V  ■  ( GVp)}dV  = 


< v 


S  L 


P 


dG 

d'Xr, 


G 


dp 

dx„ 


dS  =  0 


where  xn  is  the  vector  outward  normal  to  the  surface  S.  The  acoustic  pressure 
p  is  then  solved  for  by  substitution  into  the  integrand  of  the  Green’s  function  G 
and  the  boundary  conditions  on  S  specific  to  the  problem’s  geometry  and  source 
parameters. 


Appendix  C 


Equipment  specifications 

This  appendix  summarizes  the  important  characteristics  of  the  sensors  used  in  this 
study. 


C.l  Flow  rate  sensor,  Sierra  824  Top-Trak  Mass  Flow 
Meter 

Mass  flow  rate  range:  0  to  50  slprn 

Power  requirements:  12-15  VDC  nominal,  100  mA  maximum 
Accuracy:  ±  1.5%  of  full  scale 
Output  signal: 0-5  VDC 

C.2  Pressure  sensor,  Omega  PX209  Pressure  Transducer 

Pressure  range:  0-60  psi 
Power  requirements:  12-15  VDC 
Accuracy:  ±0.25% 

Output  signal:  0-5  VDC 
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C.3  Temperature  sensor,  Omega  SCASS-1256-6  Thermo¬ 
couple 

Temperature  range:0  —  220°  C 
Power  requirements:  12-15  VDC 
Accuracy:  ±1°  C 
Output  signal:  0-5  VDC 


C.4  In  air  microphone,  Bruel  and  Kjaer  4938  Series 
Microphone 


Dynamic  range:  30-172  dB 
Sensitiviy:  16  mV/Pa 


C.5  Unsteady  force  transducers,  PCB  Piezotronics  208C01 
and  208 All  Dynamic  Force  Sensors 

Force  input  range:  0-10  lbs 

Power  requirements:  20-30  VDC  with  constant  current  of  2-20  mA 
Sensitivity:  500  mV/lb 
Accuracy:  ±  1  % 

Output  signal:  0-5  VDC 

Internal  resonance  frequency:  32  kHz 
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C.6  Hydrophones,  G.R.A.S.  10CC  Hydrophone 


Frequency  range:  5Hz-120kHz 
Sensitivity:  -170  dB  re  1  V /mu Pa 
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